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ABSTRACT 


When  investigating  the  detection  performance  of  a  passive 
homing  torpedo  used  against  shallow  draft  surface  ships, 
certain  environmental  factors  such  as  the  rough  sea  surface 
and  the  bubble-dominated  inhomogeneous  layer  near  the  sea 
surface  have  to  be  considered.  This  thesis  attempts  to  gain 
some  insight  into  the  behavior  of  a  homing  torpedo  system 
during  its  critical  attack  phase,  as  well  as  getting  some 
indications  of  the  relative  importance  of  the  scattering 
mechanisms  and  the  induced  tactical  limitations.  An  ideal¬ 
ized  propagation  model  was  used  as  reference  of  comparison. 
For  a  given  sea  state  and  target  speed  the  results  stress 
the  importance  of  low  operating  frequency  as  well  as  a  high 
maximum  turn  rate.  They  also  point  to  the  importance  of 
having  a  search  depth  below  the  bubble-dominated  subsurface 
layer,  and  a  variable  speed  capability  during  the  torpedo's 
attack  phase. 
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I .  INTRODUCTION 


The  following  analysis  examines  several  factors  that 
limit  the  detection  performance  of  a  passive  homing  torpedo 
with  the  mission  objective  of  countering  shallow-draft 
targets  in  Norwegian  coastal  waters.  Generally,  these 
factors  can  be  divided  into  three  main  groups: 

-Environmental  factors  in  the  ocean 

-Electrical  mechanical  and  hydrodynamical  factors  in 
the  torpedo  system 

-The  users  tactical  situation 

The  factors  that  are  generated  in  the  ocean  itself  are  the 
principal  subject  of  this  analysis.  As  our  interest  is 
confined  to  the  layer  immediately  below  the  surface  of  the 
ocean,  later  called  the  subsurface  ocean  layer,  the  main 
factors  affecting  the  sound  propagation  are: 

-Scattering  and  absorption  due  to  the  bubble-dominated 
inhomogeneous  subsurface  ocean  layer. 

-Scattering  from  the  rough  sea  surface 
The  concentration  of  air  bubbles  and  the  roughness  of  the 
sea  surface  are  determined  by  the  windspeed.  The  effective¬ 
ness  of  both  these  scattering  mechanisms  depends  on  the 
frequency  of  the  incident  wave  and  the  geometry  of  the 
source  and  receiver.  The  following  analysis  is  limited  to 
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high  frequencies  in  the  region  of  30-60  kHz  which  are 
characteristic  of  existing  torpedo  systems.  At  high  fre¬ 
quency  and  low  grazing  angles  for  the  incident  and 
received  signals,  the  phenomena  of  "shadowing"  of  the 
surface  by  other  parts  of  the  boundary  occurs.  Under 
these  conditions,  the  effect  from  the  inhomogeneous  sub¬ 
surface  ocean  layer  becomes  increasingly  important.  Ob¬ 
viously  both  the  above  mentioned  scattering  mechanisms  are 
present  simultaneously.  Often  these  two  effects  cannot  be 
resolved  either  theoretically  or  experimentally,  as  any 
signal  with  a  finite  duration  will  be  scattered  from  the 
space  near  the  surface  simultaneously  with  that  from  the 
sea  surface  itself. 

In  order  to  adequately  describe  the  scattering  mech¬ 
anisms,  this  analysis  starts  with  a  presentation  of  the 
oceanographic  background  material  for: 

-Typical  windspeed  and  wave  height 

-Typical  sound  speed  profiles 

-Density  and  distribution  of  air  bubbles  in  the  sub¬ 
surface  layer 

-Statistical  description  of  the  sea  surface 
The  analysis  proceeds  by  separately  estimating  the  effect 
of: 

-Scattering  from  a  randomly  rough  surface 

-Scattering  and  absorption  caused  by  am  inhomogeneous 
subsurface  layer. 
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and  comparing  their  relative  importance.  The  method 
employed  for  these  estimations  is  an  approximation  that 
is  a  combination  of  both  ray  and  wave  theories.  Ray 
methods  are  used  to  follow  the  acoustic  signal  from  its 
source  to  the  vicinity  of  the  scatterer.  Wave  theory  is 
used  to  calculate  the  actual  scattering  process.  Finally, 
ray  theory  is  again  used  to  follow  the  scattered  signal 
to  the  receiver. 

An  idealized  propagation  model  consisting  of  an  iso¬ 
tropic  stratified  medium  will  be  used  as  reference  of 
comparison.  This  model  is  founded  on: 

-A  noise  source  from  a  cavitating  propeller  blade. 

-The  operational  characteristics  for  a  square  law 
detector  (ROC-curves) . 

-A  transmission  loss  model  based  on  geometrical 
spreading  and  absorption  losses  in  homogeneous  sea 
water. 

In  this  analysis,  the  passive  sonar  equation  is  used  to 
predict  the  performance  of  the  homing  system.  The  detec¬ 
tion  range  encountering  the  two  scattering  effects  will 
be  obtained  from  the  sonar  equation  and  compared  to  the 
detection  range  based  on  the  reference  model.  Thus,  the 
difference  in  ranges  at  which  the  homing  device  just 
acquires  the  target  with  and  without  scattering  is  a 
measure  of  efficiency. 
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II.  SCENARIO  AND  TORPEDO  RUN  GEOMETRY 

In  the  Norwegian  coastal  waters,  the  primary  mission  of 
a  torpedo  system  is  to  counter  an  amphibious  force  consist¬ 
ing  of  escorts,  supply  ships,  id  shallow-draft  landing 
crafts.  Typical  characteristics  of  these  three  ship 
types  are  as  follows: 

Supply  ships: 

-Displacement  5000  tons 
-Length  100  m. 

-Draft  6  m. 

-Speed  15-20  kts. 

Escort  ships: 

-Displacement  2000-3000  tons 
-Length  85  m. 

-Draft  3  m. 

-Speed  35  kts. 

Landing  craft: 

-Displacement  1000  tons 
-Length  80  m. 

-Draft  2  m. 

-Speed  18  kts. 

In  order  to  simplify  this  analysis,  moderate  sea  states  (SS  3) 
are  assumed.  Since  "moderate"  wave  heights  of  two  meters 
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are  appreciable  when  compared  with  the  two-meters  draft 
of  the  landing  craft,  the  possibility  of  an  acoustic  torpedo 
impacting  the  target  at  a  depth  of  two  meters  is  very  remote 
without  the  use  of  an  influence  exploder. 

Two  relevant  search  and  attack  schemes  will  be  con¬ 
sidered.  These  are  illustrated  in  Fig.  1  together  with  the 
operation  of  the  influence  exploder.  From  target  valida¬ 
tion  to  completion  of  terminal  attack,  the  torpedo  contin¬ 
uously  tracks  in  the  azimuth  plane.  In  Case  A,  ascent  is 
inhibited  after  enable.  For  Case  B,  ascent  is  inhibited 
after  the  torpedo  reaches  terminal  attack  depth. 

An  assumed  attack  depth  of  six  meters  is  consistent 
with  the  activation  range  of  influence  exploders  and  is 
deep  enough  to  preclude  wave  or  "free  surface"  induced 
disturbances  of  the  torpedo.  Success  of  the  attack  depends 
primarily  on  the  availability  of  maintaining  azimuth-plane 
steering  to  within  a  short  horizontal  range  of  the  target, 
and  the  subsequent  operation  of  the  influence  exploder. 

Case  A  is  of  particular  interest  to  this  analysis,  as 
both  the  search  and  attack-depth  are  within  the  subsurface 
layer. 


III.  OCEANOGRAPHIC  BACKGROUND  MATERIAL 


The  oceanographic  background  for  predicting  typical 
and  extreme  conditions  of 
-wind  speed 
-wave  heights 

-bubble  densities  and  distributions 
-ambient  noise  versus  self  noise 
-sound  speed  profiles 

are  outlined  in  detail  in  Appendix  A.  Even  though  most 
data  have  beneral  validity  for  Norwegian  coastal  waters, 
the  region  above  68°N  are  of  particular  interest.  Thus, 
a  typical  area  combining  open  as  well  as  confined  waters 
can  be  represented  by  "Andfjord"  at  70°N,  where  the  oceano¬ 
graphic  conditions  can  be  related  to  the  weather  station 
"Andenes,"  see  Pig.  2. 

Figures  3  and  4  [Ref.  1]  show  average  windspeed  and 
the  occurrence  of  significant  wave  heights  as  a  function 
of  time  of  year  at  weather  station  "Andenes,"  respectively. 
The  bulk  of  data  is  centered  around  a  windspeed  of 
Beaufort:4-5  (11-21  kts)  and  SS:3-4  (significant  wave  heights 
1-2  m)  .  Table  I  gives  the  relationship  between  SS ,  wind- 
speed  and  expected  significant  wave  heights. 
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TABLE  I 


RELATIONSHIP  BETWEEN  SEA  STATE  (SS),  WIND  SPEED, 

AND  EXPECTED  SIGNIFICANT  WAVE  HEIGHTS  (H  ) 

s 


Class 

Number 


Significant  Wave  Height 
in  m 


Beaufort 

Scale 


SS 


H 


s 


0 

0 

0 

1 

0  -  0.1 

1 

2 

0.1  -  0.5 

2 

3 

0.5  -  1.25 

3-4 

4 

1.25  -  2.5 

5 

5 

2.5  -  4.0 

6-7 

6 

4.0  -  6.0 

8 

7 

6.0  -  9.0 

9-10 

8 

9.0  -  14.0 

11 

9 

>14 

12 
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This,  together  with  the  low  probability  of  having  an 
amnibious  operation  occurring  in  high  sea  states  (SS>5) 
justifies  the  assumption  of  moderate  sea  state  with  wind- 
speed  in  the  region  of  12  kts  and  wave  heights  of  2  m. 

The  bubble  data  distribution  taken  from  Ref.  2  was 
obtained  in  the  area  "Tromso"  -  "Bjomoya"  during  the 
period  J tine -November  1978.  These  data  correlate  very  well 
with  a  larger  body  of  data  obtained  by  H.  Medwin  [Ref.  3]. 

Figures  5,  6  and  7  [Ref.  2]  show  the  density  of  resonant 
bubbles  as  a  function  of  depth  with  windspeed  as  parameter 
for  the  12, 38  and  120  kHz.  As  seen  from  these  figures, 
the  number  of  resonant  bubbles  are  an  increasing  function 
of  frequency  and  windspeed,  and  a  decreasing  function  of 
depth.  Below  a  depth  of  approximately  15  m,  the  number  of 
bubbles  is  negligible  for  the  windspeed  of  interest. 

The  effect  of  SS  (windforce)  on  the  ambient  noise  level 
is  given  in  Fig.  8  [Ref.  4],  Shallow  coastal  Norwegian 
waters  are  typically  5-10  dB  noisier  than  the  corresponding 
deep  water.  However,  great  variability  caused  by  local 
ship  traffic,  fishing  fleet  activity,  marine  life  and 
local  wind  conditions  makes  ambient  noise  prediction  diffi¬ 
cult  in  these  areas.  This  means  that  for  accurate  modeling, 
ambient  noise  prediction  have  to  be  done  at  each  location 
as  its  level  is  both  site  and  time  dependent.  However, 

Fig.  8  shows  that  for  frequencies  higher  than  50  kHz,  the 
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effect  of  windforce  on  the  ambient  noise  level  decreases 


to  a  lower  bound  determined  by  the  thermal  agitation. 

Based  on  the  above  discussion  and  experience  related  to 
noise  levels  for  torpedo  systems,  the  self  noise  will  be 
assumed  to  be  dominant  through  this  analysis. 

Figure  9,  obtained  from  Ref.  5,  shows  that  the  sound 
speed  profiles  usually  encountered  in  the  area  of  interest 
results  in  extremely  difficult  sonar  conditions.  This  is 
illustrated  in  Figs.  10  and  11,  which  show  worst-case  ray 
path  derived  from  Fig.  9.  In  addition,  the  presence  of 
bubbles  in  the  subsurface  layer  causes  the  sound  speed  to 
be  a  function  of  frequency.  The  above  two  factors  may 
frequently  be  the  ones  limiting  the  detection  range  of  the 
torpedo.  These  effects  can  be  minimized  by  selecting  an 
appropriate  search  depth  for  a  particular  sound  speed 
profile.  In  addition,  for  Case  B  the  corresponding  curved 
homing  trajectories  in  the  pitch  plane  give  an  error  in 
apparent  range  to  target.  This  effect  will  normally  be 
taken  into  account  by  devising  appropriate  attack  logic 
which  is  outside  the  scope  of  the  present  analysis. 
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IV.  THE  PASSIVE  SONAR  EQUATION 


A  measure  of  efficiency  for  a  passive  homing  torpedo 
is  the  detection  range  obtained  by  solving  the  passive 
sonar  equation  for  broadband  noise: 

SL+lOlogB-TL  (geom)  -aR+DI-NL-DT=*0  (1) 

where 

SL  «  spectral  level  of  the  broadband  noise 

radiated  by  the  target  (in  dB  re  luPa/lHz  at 
lm)  . 

R  =  detection  range  (in  m) . 

a  =  attenuation  coefficient  at  the  center  frequency 
(in  dB/m) . 

DI  =  receiving  sensitivity  (directivity  index) 

(in  dB  re  lyPa) . 

NL  ■  noise  level  at  the  receiver  in  the  bandwidth  B 
(in  dB  re  lyPa) . 

DT  =  detection  threshold;  the  signal  to  noise  ratio 
at  the  transducer  output  required  for  a  de¬ 
tection  probability  of  pD  and  associated  false 
alarm  probability  pFA  (in  dB) . 
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V .  REFERENCE  MODEL 


A.  INTRODUCTION 

In  order  to  produce  the  reference  for  the  analysis 
the  sonar  equation  is  solved  assuming  ideal  free-field 
conditions,  a  simple  noise  source  model,  and  a  generalized 
square-law  detector. 

B.  IDEALIZER  TRANSMISSION  LOSS  MODEL 

Because  the  presence  of  refraction,  scattering,  and  of 
ocean  boundaries,  free-field  conditions  associated  with 
homogeneous  (isovelocity)  and  unbounded  medium  seldom 
exist  in  the  sea.  However,  as  a  basis  for  comparison, 
the  ubigutuous  spherical  spreading  law  plus  an  added  loss 
term  due  to  "normal  absorption"  can  be  used  as  a  reference 
model  for  measuring  the  effects  of  the  previous  mentioned 
scattering  and  absorption  mechanisms.  Thus,  the  reference 
transmission  loss  model  can  be  expressed  as: 

TL=201ogR  +  aR  (2) 

where  the  absorption  coefficient,  expressed  in  dB/m,  can 
be  obtained  from  Fig.  14  taken  from  Ref.  8. 
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C.  NOISE  SOURCE  MODEL 

1 .  General  Characteristics  of  Noise  Sources 

Sound  is  generated  in  a  fluid  medium  by  any  process 
that  causes  an  unsteady  pressure  field.  Physically  processes 
that  can  cause  an  unsteady  pressure  field  include: 

-Pulsation  of  a  boundary  surface  of  the  medium 
-The  action  of  a  nonsteady  source  on  the  fluid 
-Turbulent  motion  in  the  fluid 
-Oscillatory  temperatures 

It  can  be  shown,  e.g..  Ref.  10,  that  each  source  mechanism 
mathematically  corresponds  to  a  dominant  order  of  multipole. 
If  all  sources  are  of  such  a  nature  that  their  time  variation 
can  be  described  by  a  Fourier  Integral,  it  can  be  shown 
[Ref.  10]  that  the  Helmholtz  Equation  is 

V2Poj(x)  +  gy  pw(x)  =  -4irfw(x)  (3) 


3Q(x) 

CO _ 

at 


+  7*F  (x)  - 

a) 


Term  1  Term  2 

where  the  right  hand  side  describes  distributed 
The  terms  on  the  right  hand  side  of  Eq. 
following  interpretations : 

Term  Is  mass  injection 
Term  2:  external  force 
Term  3:  turbulent  shear  stress 


a2T 


0) 


±1 


3xi3Xj 


Term  3 

source  terms. 
( 3 )  have  the 
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In  the  long  distance  and  long  wavelength  approximation, 
it  can  be  shown  that  the  mass  injection  term  gives  rise  to 
a  simple  source;  a  zero  order  pole  called  a  monopole.  The 
monopole  radiates  omnidirectional  with  no  angular  dependence. 
At  large  distances  the  pressure  field  from  the  monopole 
radition  is  that  of  a  point  source.  Examples  of  this  are: 
-Pulsating  bubbles 
-Cavitation 

The  external  force,  in  the  long  distance  and  long  wave 
length  approximation,  is  associated  with  a  dominant  dipole 
which  has  a  cosine  directional  pattern.  Examples  of  this 
type  of  radiation  is  that  caused  by  the  vibratory  motion 
of  an  unbaffled  rigid  body. 

Radiation  from  turbulent  shear  stresses  is  charac¬ 
terized  by  a  lowest  order  term  of  quadrupole  nature. 

The  efficiency  of  the  source  terms  decreases  with 
increasing  dependence  on  the  spatial  derivatives.  This 
can  be  understood  when  recognizing  that  wave  functions  of 
the  general  form  f(x-ct)  have  a  time  derivative 

|-~f(x-ct)|  =  cffc(x-ct)  (4) 

which  is  magnitude  c  (sound  speed)  greater  than  the  spatial 
derivative 

|^~ f(x-ct)j  =  fx(x-ct)  (5) 
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Other  factors  being  equal,  the  radiation  from  an  external 
force  is  small  compared  to  that  from  mass  injection,  and 
that  from  turbulent  shear  stress  is  the  smallest;  therefore, 
monopole  radiation  is  the  dominant  term. 

Propeller  cavitation,  when  it  occurs,  is  usually 
the  dominant  noise  source  for  any  marine  vessel.  Submarine 
and  torpedoes  often  operate  at  a  depth  great  enough  to  avoid 
cavitation.  Surface  ships,  on  the  other  hand,  generally 
have  well  developed  propeller  cavitation  with  the  result 
that  their  radiated  spectrum  from  5  Hz  to  100  kHz  is 

controlled  by  this  source. 

The  basic  phenomena  of  cavitation  combined  with 

propeller  hydrodynamics  give  the  fundamental  characteristics 
of  propeller  cavitation  noise.  An  excellent  qualitatively 
discussion  of  this  can  be  found  in  Ref.  9:Chs.  7  and  8, 
from  which  the  following  is  extracted: 

Propeller  blades  are  rotating  twisted  wings 
that  produce  hydrodynamic  forces.  Depending  on  operating 
conditions,  they  experience  cavitation  on  a  number  of 
different  places.  Of  these  there  are  three  prominent  types: 
-Tip  vortex  cavitation 
-Hub  vortex  cavitation 
-Surface  blade  cavitation 
In  addition  to  the  two  types  of  vortex  cavitation, 
there  normally  are  two  types  of  blade  surface  cavitation: 
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-Back:  driving  face 
-Front:  suction  surface 

Of  all  kinds  of  propeller  cavitation,  surface  blade  cavi¬ 
tation  on  the  suction  surface  is  normally  the  most  noisy, 
while  hub  vortex  cavitation  is  the  least  noisy. 

2 .  The  Noise  Source  Model 

Due  to  lack  of  recorded  and  available  noise  data 


from  the  target  in  question,  the  noise  source  has  been 
generalized  on  the  basis  of  the  following  discussion  and 
assumptions: 

The  noise  source  will  be  build  up  around  a  surface 
blade  cavitating  propeller  operating  in  a  good  to  poor 
wake;  surface  cavitation  will  be  assumed  to  be  dominant. 

D.  Ross  [Ref.  9]  has  developed  an  approximate 
theory  for  cavitation  noise,  where  dimensional  analysis 
is  combined  with  the  basic  results  of  cavitation  theory 
that  the  acoustic  pressure  is  proportional  to  the  product 
of  the  collapse  pressure  of  the  cavities  and  the  volume 
of  cavitation  produced  per  unit  time.  From  this  synthesis 
it  is  found  that  the  total  acoustic  intensity  varies  as: 


I  *  K.. 
ti 


— t'irrXu11-  -1*'! 


where 


r  =  distance  of  the  hydrophone  from  the  source 
K. .  =  the  cavitation  inception  parameter 
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Ut  =  blade  tip  speed 

Uti  =  blade  tip  speed  for  inception  of  cavitation 

This  expression  shows  that  propeller  cavitation  noise  power 
is  proportional  to  the  total  number  of  blades,  b,  the  blade 
chord,  s,  and  to  the  propeller  diameter,  D,  and  is  a  function 
of  the  tip  speed  with  the  dependence  on  the  tip  speed  being 
the  strongest.  The  different  blade  surface  sections  where 
cavitation  exist  are  uncorrelated  and  the  radiated  noise 
is  treated  as  a  single  monopole  radiation  so  that  at  a 
distance  r>>a  (where  a  is  the  characteristic  dimension  on 
the  source  region)  the  radiation  is  similar  to  that  of  a 
point  source  with  no  angular  dependence. 

Submarines  and  torpedoes  with  centerline  propellers 
have  a  relatively  symmetric  inflow  condition.  Surface  ship 
propellers,  in  contrast,  operate  under  nonuniform  inflow 
conditions.  Circumferential  wake  variation  causes  the 
radiated  sound  to  be  amplitude  modulated  at  the  blade  rate 
frequency.  Furthermore,  slight  physical  difference 
between  the  blades  produces  modulation  at  the  shaft  rate 
frequency.  These  amplitude  variations  gives  a  very 
distinct  characteristic  to  the  radiated  noise  that  can  be 
used  for  classification  purposes  to  reduce  the  probability 
of  false  alarm. 

The  most  complete  source  of  data  on  surface  ship 
radiated  noise  are  measurements  made  during  WWII,  reported 
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in  a  compendium  issued  by  the  U.S.  Office  of  Scientific 
Research  and  Development  (OSRD)  in  1945  and  declassified 
in  1960.  When  these  data  are  examined  the  radiated  noise 
is  found  to  depend  on  tip  speed  and  the  number  of  propeller 
blades  with  little  dependence  on  the  other  variables.  For 
surface  ships  near  cruise  speed,  the  source  level  for 
frequency  over  100  Hz  can  be  written  as: 

SL=SL ' +2 0-20 log f ;  f>100  Hz  (7) 

where 

f  =  frequency  in  Hz 

SL' =  overall  source  level  in  dB  re  luPa. 

The  overall  level  can  be  expressed  as: 

SL'=175+601ogU  /25  +101ogb/4  (8) 

where 

U  =  rrnD  (9) 

n  =  rotational  speed  (rpm) 

D  »  diameter  of  the  propeller  (m) . 

The  above  expressions  are  used  as  the  basis  for  the  noise 
model  with  the  following  input  data: 

n  **  300  rpm  for  maximum  cruise  speed  of  15  kts. 
n  =  180  rpm  for  a  cruise  speed  of  10  kts. 

D  *  2  m. 
b  *  5. 

The  resulting  noise  spectrum,  in  dB  re  lyPa  at  1  m, 
as  a  function  of  speed  in  kts.  are  tabulated  in  Table  IX. 
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SOURCE  LEVEL  DATA  AND  CALCULATION 
A  BLADE  SURFACE  CAVITATING  PROPELLER 
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and  plotted  in  Fig.  15.  A  one  sigma  region  (5  dB  of  un¬ 
certainty)  is  incorporated  in  Fig.  15. 

The  above  noise  spectrum  estimation  agrees  very 
well  in  the  high  frequency  limit,  with  more  recent  studies 
by  A.  Lovik  [Refs.  11  and  12].  Here  it  is  found  that  the 
cavitation  spectra,  both  theoretically  and  experimentally, 
can  be  divided  into  four  frequency  regions,  as  illustrated 
in  Fig.  16. 

Region  1  is  dominated  by  noise  at  the  blade  fre¬ 
quency  and  its  harmonics.  The  emitted  sound  is  caused  by 
the  volume  variation  of  the  main  cavity. 

Region  II  starts  at  the  bubble  frequency,  which  is 
the  reciprocal  of  the  lifetime  of  the  main  cavity.  The 

mean  power  level  is  found  to  decrease  with  increasing  fre- 
*—2.5 

quency  as  f 

Region  III  is  am  inter  .  ; : ate  region. 

Region  IV  associated  with  the  shock  waves  starts 

at  the  mean  collapse  frequency  fQ,  given  by  the  mean 

collapse  time.  The  power  level  is  found  to  decrease  as 
_2 

f  ,  as  in  Eq.  5. 

The  number  of  gas  bubbles  in  the  water  have  a 
pronounced  effect  on  the  high  frequency  cavitation  noise 
from  the  propeller.  This  is  illustrated  in  Fig.  17, 
obtained  from  Ref.  12,  where  the  power  is  found  to  decrease 
as  much  as  40  dB  with  increasing  gas  content. 
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Scaling  laws  are  developed  [Ref.  12]  for  each  region 
based  on  a  series  of  models  and  full  scale  measurements. 

These  laws  depend  on  the  dynamic  pressure  induced  by  the 
propeller,  the  model  ratio,  and  the  gas  content  of  the  water. 
The  full  scale  measurements  were  performed  in  cooperation 
with  the  Royal  Norwegian  Navy  and  the  Marine  Institute  of 
Norway.  The  model  experiments  were  performed  in  the  largest 
cavitation  tunnel  at  the  Ship  Research  Institute  of  Norway. 

In  summary,  the  scaling  of  cavitation  noise  was 
demonstrated  to  be  a  useful  tool  in  predicting  a  full  scale 
cavitation  noise  as  shown  in  Fig.  18  [Ref.  12]  which  compares 
measured  noise  spectra  for  the  model  and  full  scale 
measurements. 

For  the  high  frequency  region,  the  source  levels 
are  of  the  same  magnitude  as  predicted  by  the  WWII  empirical 
formula. 

D.  PASSIVE  MODE  RECEIVER  CHARACTERISTICS 
1.  Assumptions 

For  receiver  characteristics  assume  a  square  law 
detector  with  a  center  frequency  f-60  or  30  kHz  and  a  band¬ 
width  B.  The  detection  scheme  is  shown  in  Fig.  16.  The 
principal  assumptions  employed  in  the  derivation  are  as 
follows: 

-Gaussian  signals  in  Gaussian  noise 

-Frequency  independent  signal  and  noise  spectra 
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-Integration  time  T  is  sufficiently  long  to 
permit  application  of  the  central  limit  theorem. 

2 .  Derivations 

The  detector  input  r(t)  is  assumed  to  be  a  zero- 
mean  Gaussian  process  composed  of  noise  alone  or  signal 
plus  noise  expressed  by  the  two  well  known  hypotheses 

f^rtt)  =  n  (t)  (10) 

Hj  r  (t)  =  s  (t)  +n  (t) 

where : 

n(t):  noise  signal 
s(t):  signal. 

The  two  signals  s (t)  and  n(t)  are  assumed  to  be  independent. 

Assume  that  the  spectral  shape  of  s(t)  and  n(t) 
are  the  same,  such  that  Hq  and  only  differ  in  the  total 
power  level.  Then  the  detector- smoother  have  the  form: 


VJ 

Out 

J 

In 

T 

ifxUUi  _ 1 


Schematic  of  detector-smoother 

and 

x (t)  *  r2 (t)  (11) 

Furthermore,  let  the  noise  variance  be  normalized  to  unity 
(for  convenience)  and  the  signal  variance  by  denoted  by  a2 
Var[n(t) ]«1  (12) 

Var [s (t) ]«o2 
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Because  of  the  assumed  similarity  in  the  spectral  shapes, 
the  autocorrelation  functions  are 

Rn  ( t )  =F  ’[N  (f )  ]  =  p  (t )  ;  N  ( f )  is  the  noise  power  (13) 
spectral  density. 

Rg  (  t  )  =  p  (  t  )  a  2 
Rn(t)=  Hq  p(x) 

H]_  (l+o2)p(x) 

Furthermore,  assume  that  the  integration  time  T  is  long 

enough  so  the  central  limit  theorem  holds,  implying  that 

K  also  is  a  Gaussian  random  variable. 

This  yields  that  the  probability  density  function 

of  the  output  variable  and  hence  the  detection  and  false 

alarm  probabilities  are  completely  determined  once  the 

mean  and  the  variance  of  K  are  derived. 

If  a  process  V(t)  is  wide-sense  stationary,  then 
T  T 

E[V(t)]=E|  fv(t)dt  =|  ^/e [V ( t)  ]  =  v  (constant)  (14) 
0  0 


Thus,  assuming  that  r(t)  is  a  wide-sense  stationary  process. 

E[K]=E[x(t) ]=E[r (t) ]=l+a2  (15) 

and  similarly 


Var [V]=E [V2 ] -(E [V] }  2=  -k- 

T 


TT 

jj r  E[V(t)V(s)  ]dtds— v2 
00 


(16) 


TT  TT 

Var[V]=*-i-  jj  [Ry(t-s) -v2  ]dtds*^f  jj  Covv  (t-s)  dtds 
00  00 
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where  Cov^(t-s)  is  a  covariance  function. 
Then  letting 

T  =  t-S 


z,  =  t+s,  and  substituting  into  Eq.  (7)  yields 

15 


rr. 


Var  t  V1  =^r// Covv ( T  > 


r  t 


Var(V>=i  /  [l-i^]Covv(T)dT 


-T 


Consequently 

T 

Var[K]=i  f  [l-M]Covv  (t)  dr 
-T 

Then,  evaluating  the  covariance  function  from  the  auto 
correlation  function 

Rx(T)=E[x(t)x(t-x) ]  =  E [r2 ( t) r2 (t-x ) ] 

Since  r(t)  is  Gaussian,  the  above  fourth  moment  can  be 
expressed  as  product  and  sums  of  second  moments: 


Rx(T)=»r2(t)  *rz(t+T)  +  2r(t)r(t-x)  *r(t)r(t-x) 

=R2(0)  +  2R2(x) 

Rx(x)=(l+a2) 2+2 (l+o2)2  p  2  ( t ) 

Thus,  the  covariance  function  is 

Covx(T)“Rx(T)"{E[x(t) ^  ^ 

»(l+a2)  2+2(l+a2} p2  (t)-(H-o2)  2 


(17) 


(18) 


(19) 


(20) 
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Cov^ ( t ) =2 ( 1+a  )  p 2 ( t ) 

Inserting  Eq.  (12)  into  Eq.  (9)  yields 


X 

Var[K]=-|-  /  { [i— Li-L]  2  (1+a2 )  2 p2  (x)  }  dx 
-T 

Var[K]=iii^il  f  [1-M]  P2<^ 


If  T  is  large  compared  to  the  correlation  time 
TB>>1 

then  we  can  substitute  the  limit  for  Eq.  (22)  by 


Var  [K]  = 


2  (l+o2 ) 


p  2 ( t )  dx 


2 (1+a2 ) 2 
T 


N2  (f)df 


If  we  further  make  the  assumption  that  the  signal  and  noise 
have  ideal  flat  bandpass  spectra: 

N  ( f )  =  1/2B ,  f-B/2<f <f+b/2  (24) 

k 0 ,  otherwise , 

Inserting  Eq.  (24)  in  Eq.  (23)  yields 


Var  [K]  = 


(1+a2) 


The  probability  density  function  for  the  output  of  the 
detector  have  the  following  form 
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Here  w  is  the  outcome  of  all  possible  signals.  The  false 
alarm  probability  is  obtained  by  integrating  the  conditional 
probability  p(w|HQ)  over  the  outcome  space  for  which  to 
choose  . 


PFA 


l  p(w|HQ)dw  =  0(^> 


(26) 


Similarity,  the  miss  probability^l-detection  probability 
is  obtained  by  integrating  the  conditional  probability 
p(w|H^)  over  the  outcome  space  for  which  to  choose  Hq. 


•M 


*  1_PD  = 


.  U,  "j 

p(w|H1)dw  =  Q  ( — — — ) 


a  i 


(27) 


Further  defining  the  input  and  output  signal-to-noise 
ratios  as 

2 

S/N (input)  =  =  o2  (28) 


(Ui-U,) 2 

S/N  (output)  =  - - -  (29) 

a  i 
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For  the  square  law  detector  where 


UQ  =  E [n2 (t) ]  =1,  =  E [r2 (t) ]  =E[K]  =l+o2  (30) 

° 0 2  =  Var [n2  (t)  ]  =  Ql2  =  Var[r2(t)] 


=  VartK]  - 


The  output  signal-to-noise  ratio: 


S/N  (output)  =  -1--  ■  =  5^-2 -  (31) 

rjl+oV,  (1+a2) 2 
L  BT  J 

The  input  signal-to-noise  ratio: 

S/N (input)  *  a2  (32) 

The  probability  of  false  alarm: 

PFA  =  Q[/BT(j-l)  ]  (33) 

The  probability  of  detection: 

p  =  1-Q[/Bf  (— °— p-)  ]  (34) 

1+a 

The  Equations  (32),  (33),  and  (34)  are  plotted  for  a 
variety  of  realistically  encountered  conditions.  Figure  17 
gives  BT  versus  S/N  (in)  (identical  to  the  detection 
threshold  DT)  for  various  combinations  of  pD  and  PFA* 

Figure  18  gives  pQ  versus  S/N (in)  for  various 
combinations  of  the  threshold  j  and  PFA*  This  constitutes 
the  ROC -curves  for  the  generalized  square  law  detector. 
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E.  REFERENCE  DETECTION  PERFORMANCE 

1 .  Introduction 

As  both  the  scattering  and  absorption  are  frequency 
dependent,  it  is  necessary  to  have  reference  models  for  both 
60  and  30  kHz. 

2 .  60  kHz  Case 

The  range  dependent  portion  of  the  sonar  Equation  (1) 

-201ogr-ctR 

is  plotted  in  Fig.  19  for  a=0.021  dB/m  taken  from  Fig.  14. 

Figures  15,  17  and  19  are  then  used  to  estimate  the 
detection  range : 

(a)  Select  pQ,  PFA»  B,  and  the  integration  time  T. 
The  detection  threshold  DT=S/N (input)  is  then 
found  from  Fig.  17. 

(b)  Select  the  speed  of  the  target  and  find  the  SL 
from  the  noise  source  model  (Fig.  15) .  Then, 
reasonable  values  for  the  receiver  sensitivity 
DI  and  the  self  noise  level  NL  yields  the  left 
hand  side  of  the  sonar  equation  (1)  except  for 
range  dependent  term. 

(c)  Use  Fig.  19  to  solve  the  passive  equation 
for  R. 

A  realistic  example  may  illustrate  the  above  procedure. 
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} 


(a)  Entering  Fig.  17  with 


PD  =  0.5 


B  =  4500  Hz 
T  =  100  msec 


yields 

DT  =  -6.5  dB. 

(b)  Selecting  a  target  speed  of  12  kts.  gives 
(from  Fig.  15)  SL=100  dB.  Selecting  a 
typical  transducer  sensitivity  DI  =  -180  dB. 
Assuming  the  NL  to  be  dominated  by  self  noise 
of  typical  value  NL=-124  dB.  This  yields 

SL+101ogB+DI-NL-DT=87  dB. 

(c)  Figure  19  with 

-201ogR-0 . 02R=-87 
yields 

R=1200  m.  for  a=0.02  dB/m  and  f-50  kHz. 

The  influence  of  different  design  parameters  like  self 
noise  and  detection  threshold  on  the  passive  detection 
performance  is  now  easily  investigated  by  the  above 
procedure . 

Although  outside  the  main  scope  of  this  analysis, 
the  above  statement  can  be  confirmed  with  an  example.  As 
seen  from  Fig.  18,  a  probability  of  false  alarm  ppA  =  10  6 
implies  a  threshold  setting  =  0.9  dB  above  the  noise  level. 
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Due  to  the  variability  of  the  noise  level  together  with 
the  practical  difficulty  in  accurately  setting  the  threshold, 
a  more  realistic  goal  for  the  threshold  would  typically  be 
3  dB.  Going  into  Fig.  17  shows  that  the  corresponding  value 
for  DT  for  PD  =  0-5  is  DT  =  0  dB,  with  a  corresponding  low 
value  for  the  PFA- 

Letting  DT=0  dB  and  keeping  the  previous  assumed 
values  of  SL,  DI,  and  NL  yields: 

SL+101ogB+DI-NL-DT=80. 5  dB. 

The  corresponding  detection  range  is: 

R  =  1000  m,  for  a  =  0.02  dB/m  and  f  =  60  kHz. 

Thus,  this  change  in  threshold  setting  caused  a  decrease 
in  detection  range  from  1200  m  to  1000  m  in  return  of  a 
significant  decrease  in  the  false  alarm  probability. 

3 .  30  kHz  Case 

In  order  to  estimate  the  reference  detection  range 
for  an  operating  frequency  of  30  kHz,  we  utilize  the  sonar 
equation  (1) . 

Assume  that  the  receiver  has  the  same  generalized 
passive  detector  characteristics  as  in  the  60  kHz  case: 

DI  =  -180  dB 

NL  =  -124  dB. 

DT  =  -6.5  dB,  based  on  pD  =  0.5  and  p?A  =  10  6. 
However,  the  empirical  equation  (7)  for  the  source  level: 

SL*SL'+20-201ogf 
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—  2 

shows  that  the  source  level  falls  off  as  f  .  If  the 
dynamical  and  dimensional  parameters  of  the  propeller  are 
the  same,  SL  will  increase  by  +6  dB  when  the  frequency  is 
reduced  from  60  to  30  kHz. 

A  source  level  of 
SL=100+5=106  dB 

gives  a  range  dependent  solution  of  the  sonar  equation 
-201ogR  -otR=-87-6  =  -93  dB. 

A  plot  of 
-201ogR-aR 

is  given  in  Fig.  20  for  an  absorption  coefficient  a=0.01  dB/m 
taken  from  Fig.  14. 

Figure  20  then  gives  a  detection  range  of 
R  =  2400  m. 

Thus,  as  seen  from  these  ideal  reference  calculations, 
halving  the  frequency  gives  a  higher  source  level  and  a 
lower  absorption  loss,  resulting  in  a  doubling  of  the 
detection  range. 
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VI .  THE  EFFECT  OF  SURFACE  SCATTERING 


A.  OCEANOGRAPHIC  DESCRIPTION  OF  THE  SEA  SURFACE 

The  roughness  of  the  sea  surface  is  the  essence  of 
the  scattering  mechanism.  Thus,  to  adequately  describe 
the  scattering  of  sound  from  a  randomly  rough  sea  surface, 
it  is  necessary  to  formulate  a  suitable  description  of  the 
sea  surface  from  an  acoustical  propagation  point  of  view. 

Generally  the  shape  of  the  rough  sea  surface  is  most 
appropriately  described  in  terms  of  time  and  spatial 
dependent  random  variables.  However,  observation  of  the 
ocean  under  the  same  environmental  (meteorological)  condi¬ 
tions  indicates  that  the  roughness  is  the  same  over  large 
areas  and  for  periods  of  at  least  several  hours.  The 
random  processes  responsible  for  the  structure  of  the  sea 
surface,  therefore,  can  be  considered  stationary  at  least 
over  periods  of  hours.  With  this  assumption,  the  sea 
surface  can  be  described  in  terms  of  the  statistical 
description  of  the  surface  displacement  function,  and  the 
distribution  of  signals  reflected  from  the  sea  surface 
can  then  be  related  to  this  probability  distribution 
function. 

The  most  significant  statistical  parameters  describ¬ 
ing  the  scattering  mechanism  from  the  randomly  rough 
surface  are: 
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-the  mean  square  slope 

-the  mean  square  surface  height 

-the  correlation  length. 

Optical  measurements  made  at  sea  by  C.  Cox  and  W.  Munk 
[Ref.  13]  showed  that  the  sea  surface  with  an  arbitrary 
wide  continuous  spectrum  of  waves  is  characterized  by  a 
Gaussian  distributed  surface  slope.  The  mean  square 
slope,  determined  from  these  optical  measurements  is 

<Cl2>  -  z*  =  ( 3+5 . 12w) xlO” 3  (35) 

where 

W=wind  speed  in  m/s  measured  41  ft  (12.5  m)  above 
the  sea  surface. 

The  Gaussian  distribution  of  surface  slopes  implies  that 
the  surface  displacement  function  can  be  described  by  a 
Gaussian  probability  density  function  with  zero  mean 

<C>  =  0 

and  variance 

<C>  =  a2 

and  Gaussian  correlation  function 
<C  (t)  t,  (t+x)  > 

The  mean  square  height  a2,  is  obtained  by  integrating  the 
frequency  spectrum  of  the  fully  developed  sea.  The  fre¬ 
quency  spectrum  G  is  given  by  the  commonly  accepted  semi- 
empirical  expression  of  Piersom-Moskowitz  [Ref.  14]  as 
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(36) 


where 


G(ft)  —  —2—  exp  [  —  6  "t?1] 

ft 5 

ft  =  frequency  (in  s 

a  =  8.1 


8  =  .74 

ftQ=  g/W  (in  s-i)  . 

W = wind  speed  in  m/s  at  19.5m  above  the  sea  surface 
g = gravitational  acceleration  (in  m/s2). 

This  gives 

00 

a2  =  /  G (ft) dft  =  HliL-  (37) 

<f  4Sg2 


For  a  Gaussian  autocorrelation  function  expressed  as 


■Mt)  =  —  <c(tK(t+x) 
a2 


>-  r‘ 


(38) 


where 

T  =  correlation  length. 

The  following  relationship  for  the  mean  square  height  holds 
for  sea  of  small  roughness 

I2  -  —  ,  see  later  Eq.  (60)  (39) 

T  2 

and  the  correlation  length  is  thus: 


T  -  *  l 


(40) 
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B.  SCATTERING  THEORY 

All  real  boundaries  are  rough  for  radiation  with  short 
enough  wavelength,  and  the  apparent  roughness  depends  on 
the  "viewing"  conditions.  The  wave  reflected  by  a  plane 
surface  has  the  same  properties  as  the  incident  wave  since 
the  radiation  is  scattered  coherently  and  there  is  a  definite 
relation  between  the  incident  and  scattered  waves. 

A  randomly  rough  surface,  however,  such  as  the  wind 
generated  ocean  surface,  scatters  radiation  in  all  directions, 
i.e.,  an  illuminated  area  is  visible  from  any  direction. 

Heuristically  there  are  two  distinct  approaches  to  this 
phenomena . 

1.  If  the  boundary  is  rough  most  of  the  radiation  is 
scattered  and  there  is  little  transmission  in  the 
specular  direction.  Thus,  the  attenuation  caused 
by  the  irregularities  can  be  included  in  the  trans¬ 
mission  equation. 

2.  If  the  surface  is  truly  smooth,  it  can  be  assumed 
that  the  effect  of  the  boundary  is  to  supplement  the 
original  pressure  field  by  an  out-of-phase  image 
contribution.  For  a  randomly  rough  surface  the 
reflected  sound  neither  completely  cancels  the 
direct  sound  nor  adds  to  give  +6  dB  pressure  peaks 
of  the  interference  pattern.  For  a  rough  surface, 

this  supplement  is  a  small  fraction  of  the  direct 
path. 
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The  second  approach  will  be  used  with  the  simplified 
assumption  that  the  sea  below  the  surface  has  an  isotropic 
statistical  description;  i.e.,  the  mean  acoustic  velocity 
and  the  mean  density  are  assumed  to  be  constant  and  have 
negligibly  small  mean  square  fluctuations. 

The  estimation  of  the  scattering  is  based  on  an  approxi¬ 
mation  method  employing  both  ray  and  wave  theory.  Ray 
methods  are  used  to  follow  the  acoustic  signal  from  the 
noise  source  to  the  vicinity  of  the  sea  surface.  Then, 
wave  theory  is  used  to  calculate  the  scattering  process. 
Finally,  ray  theory  is  used  to  follow  the  scattered  signal 
to  the  receiver. 

The  geometry  is  given  in  Fig.  24a.  The  origin  of  the 
coordinate  system  is  at  the  center  of  the  illuminated  area. 
The  x-y  plane  coincides  with  the  mean  of  the  rough  surface 
as  averaged  over  the  illuminated  area. 

The  source  and  receiver  are  at  distances  and  R2, 
respectively,  from  the  origin.  R^  is  the  xz-plane  and 
makes  the  angle  0^  with  the  z-axis.  R2  makes  the  angle  02 
with  the  z-axis  and  the  projection  of  R2  on  the  xy-plane 
has  an  angle  9^  relative  to  the  x-axis. 

For  high  frequencies  R^  and  R2  are  much  larger  than 
the  acoustic  wavelength.  Then  both  the  incident  wave 
and  scattered  waves  can  be  treated  as  nearly  plane  waves. 
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The  formulation  of  the  scattering  problem  will  be 
based  on  the  Helmholtz  integral  which  requires  known 
values  of  the  normal  derivatives  of  the  incident  and 
reflected  waves  on  the  boundary.  These  are  estimated 
by  using  the  neuristic  Kirchhoff's  approximation,  which 
assumes  that  the  wave  is  locally  reflected  by  a  plane 
surface;  i.e.,  an  approximation  restricted  to  a  surface 
not  too  rough  and  not  shadowed. 

Further,  the  receiver  derictivity,  as  indicated  by 
Fig.  24b,  will  be  used  to  limit  the  surface  area  that  is 
illuminated. 

Since  this  procedure  is  based  on  a  detailed  develop¬ 
ment  by  I.  Tolstoy  and  C.  S.  Clay  [Ref.  15],  only  the 
main  points  will  be  outlined  here  to  bring  out  the  assump¬ 
tions  made  and  the  inherent  limitations  of  this  approach. 

The  development  starts  by  considering  the  inhomogeneous 
wave  equation  of  the  general  form; 

.  32p(xjt)  _ 

V2p(Xjt)  - - -  -4tr  f(Xjt)  (41) 

c2  at2 


where 

f(xtt)  *  is  a  known  source  distribution. 

The  development  is  based  on  the  following  initial  assumptions 
-The  medium  is  homogeneous. 

-The  medium  is  bounded  by  some  surface  S,  onto  which 
an  incident  wave  impinges. 
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-The  boundary  is  characterized  by  the  specific 
acoustic  admittance  and  the  shape  of  the  boundary. 

-The  incident  wave  is  harmonic. 

The  assumed  harmonic  source  implies  that  f(x,t)  can  be 
decomposed  into  a  Fourier  integral.  Furthermore,  assuming 
that  the  solution  of  Eq.  (41)  can  be  decomposed  in  time, 
we  arrive  at  the  Helmholtz  equation 


V2p  (x)  +  p  (x)  = 


U) 


0) 


-4?rf  (x) 

0) 


(42) 


It  should  be  noted  here  that  L.  Fortuin,  in  [Ref.  16] 
showed  that  the  Helmholtz  equation  is  not  exactly  correct 
for  a  medium  with  a  time  dependent  boundary.  The  equation 
can,  however,  be  used  with  a  good  approximation  when  the 
time  derivative  of  the  surface  elevation  is  much  smaller 
than  the  speed  of  the  waves  through  the  medium.  For  under¬ 
water  sound  waves  scattered  by  the  rough  sea  surface,  this 
means  that  the  wind  speed  has  to  be  much  less  than  the 
sound  speed;  a  requirement  easily  fulfilled  for  our  inves¬ 
tigation. 

Green's  method  allows  the  solution  of  this  linear  in¬ 
homogeneous  wave  equation  to  be  expressed  in  the 
heuristic  Helmholtz  integral  form: 


u> 


(xx)G  (xx‘)d3xl+ 

CO 


,  f  3p  (x1) 

^{G(xxM  w 


3n 1 


(43) 


3G  (xx‘) 

-  p,  (x1)— ^ - }da 1 


u> 


3n 1 


1 

1 


'i 

% 


i 

i 
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The  first  integral  on  the  RHS  of  Eq.  (43)  contains 
the  sound  sources  and  the  bulk  (volume)  scattering.  The 
second  integral  represents  the  surface  scattering  and  is 


taken  over  all  finite  surfaces. 

Now,  disregarding  the  direct  path,  the  signal  at  the 
receiver  is  given  by  the  surface  integral  alone: 


plf)(X2)=T7  AG(*2XM  — • j  -  pu)(x1)^-}da1 
>  an  dn 


(44) 


Wi^ere  the  subscript  s  denotes  the  scattered  field.  In 
order  to  solve  Eq.  (44)  the  following  must  be  done: 

-Give  an  approximate  expression  for  the  incident 
wave . 

-Find  an  appropriate  Greens  function. 

-Make  an  approximation  for  p(j(x1)  and  3p^/3n  at 
the  surface. 

As  we  already  have  assumed  a  simple  harmonic  source, 
the  incident  wave  can  be  expressed  as: 


(i) 


P 


0) 


ikR 

e 


BD  ikR 

IT  e 


(45) 


where 

II:  power  output. 

D:  illumination  function. 

Assuming  kR>>l,  i.e.,  that  the  distance  of  the  source 
is  large  compared  to  the  wavelength,  the  wave  in  the 
bounded  ensonified  area  can  be  considered  as  a  plane  wave 
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characterized  by  its  propagation  vector 


k.  =  -k- 

i  i 


thus  obtaining  the  expression  for  the  incident  wave: 

P  <i)  (I  ,  §2  eikRl  e1*1  •*'  (46) 

Further  assuming  that  the  receiver  is  sufficiently  far 
removed  from  the  scattering  area  with  the  rest  of  the  sea 
surface  not  contributing,  then  the  scattering  area  acts  as 
a  small  induced  source  in  a  free  space  and  we  can  approxi¬ 
mate  the  propagation  of  the  scattered  waves  from  the  en- 
sonified  region  in  terms  of  the  free  field  Greens  function: 


where 


G  ( X  2  X  l  ) 


k  —  k-|-=r 
s  X 


We  further  assume  that  each  surface  element  da1  acts  as  a 
small  reflector,  and  that  the  response  of  da1  to  the  inci¬ 
dent  wave  is  that  of  a  "local  reaction,"  i.e. ,  independent 
of  any  other  part  of  the  ensonified  area  A. 

Then,  the  Kirchhoff's  approximation  where  it  is  assumed 
that  p  and  3p/3n  vanish  everywhere  on  the  surface  except  at 
the  ensonified  area  and  that  the  values  of  p  and  3p/3n  are 
proportional  to  the  incident  wave,  allow  the  scattered 
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Y  =  -|(COS  0!  +  COS  02 ) 


C (x  y  )  =  surface  displacement  function. 

The  above  scattering  integral  is  then  applied  to  a  randomly 
rough  surface  where  the  surface  displacement  function  r  is 
a  random  variable  assumed  to  be  represented  by  a  Gaussian 
PDF  expressed  as 


WU)  = 


e 


-C2/2 a: 


(51) 


a/Tn 

with  zero  mean  and  variance  a2. 

Also,  assume  that  the  surface  is  slowly  varying  so  that 
the  signal  reflects  from  an  essentially  stationary  surface 
and  that  the  succession  of  received  scattered  signals  pn 
are  assumed  to  form  a  satistically  independent  set  of 
sample  functions  from  which  sequence  N  the  first  and  second 
moment  of  the  field  can  be  investigated.  Doing  so  Clay  and 
Tolstoy  [Ref.  15]  found  that  the  mean  reflected  signal  can 
be  expressed  as 

-2y 2a2 


<P, 


(s) 


(x2) >C  =  PQe 


(52) 


where 


P0  *  PoI3’  Ul) 


C=0 


is  the  signal  reflected  by  a  mirror-like  surface, 
other  factors  being  the  same. 
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It  is  seen  from  the  above  thac 

1.  For  a->-0,  all  displacements  c  have  zero  proba-  j 

bility  and  the  mean  signal  tends  to  pQ .  Further-  , 

more,  all  elements  contribute  to  the  scattering 

coherently.  j 

2.  For  a>>k,  all  displacements  £  are  equally  probable.  j 

j 

There  are  large  phase  shifts  between  contributions 

from  different  surface  elements  and  they  tend  to  cancel 

each  other  and  the  scattering  radiation  is  incoherent.  j 

The  second  moment  is  defined  as: 

T 

<s2(t)  ~  \  f  s 2 (t) dt  =  <pp*>T  -<pp*>T  (53) 

where 

1  N 

s2  (t)  =  Pn2(t)  -  N<p(t)>N2]  (54) 

In  this  expression  the  operation  of  squaring  the  signal  has 
to  be  considered.  However,  each  surface  element  has  a 
different  c  for  a  random  surface.  Thus,  the  probability 
of  finding  element  dx'  dy'  with  Z i ,  and  element  5  2  with 
dx"  dy"  is  expressed  in  the  bivariate  distribution  function 
assumed  to  be  Gaussian  and  of  the  form: 

WUiCa)- - - - r  e  [ - - -  ( ;  1  2  +  C  1  z-2c  !  C  2'H  ]  (55) 

2tto2  (1-t^2 ) **  2(1  -<J>2)o2 

where 

<M£»n)  =  —  <C  1  (x‘y 1 1)  ;2(x"y"t)>  (56) 


•i 

.* 


i 


is  the  cross-correlation  function  characterizing 
the  surface  shape. 

Changing  to  polar  coordinate  leads  to  Eq.  6.51  of  Ref.  15 
where : 

<s2 >a  /  DJ  (2Kr>  [e-4Y2a1<1-.l.)_e-4Y20s)  rdr  (57 

0  o 

Here  J  is  the  Bessel  function  of  zero  order  and  x  is  the 
o 

transformation  parameter  given  as 

K=  _  o.._ 

cos  6 

As  seen  from  Eq.  (57) 

1.  For  rough  surface,  y'2a2  is  large  and  the  second 
term  in  the  bracket,  the  coherent  part,  is 
negligible. 

2.  For  smooth  surface,  y2a2  is  zero  and  the  whole 
bracket  is  zero. 

Since  Eq.  (57)  cannot  be  integrated  directly,  Clay  and 
Tolstoy  [Ref.  15]  consider  it  for  small  and  large  yo 
separately. 

Thus,  concentrating  on  the  high  frequency  limit,  Clay 

and  Tolstoy  assumed  y2a2>>l  and  showed  that: 

-The  coherent  component  is  negligible. 

— <s 2 >  =<pp*>,  as  the  means  tends  to  zero, 

00 

<s2  >=<pp*>a  (2«r)  e_4Y  a 


rdr;  y2a2>>l  (58 
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Equation  (58)  consists  of  the  product  of  an  oscillatory 
function  and  an  exponential  function.  Because  of  the  Bessel 
function,  the  main  contribution  to  the  integral  is  near  r=0. 
Near  r=0  the  phase  changes  slowly  and  the  expression  can  be 
evaluated  by  the  method  of  stationary  phase.  Thus,  the 
expression  of  ip  about  r=0  is  given  as: 

ip  z  l+<(»"(0)^-  (59) 

Furthermore,  Clay  and  Tolstoy  show  that  ip"  (0)  can  be  related 
to  the  characteristics  of  the  surface  as: 


|  'p"  (0)  |  =  -4-  <C12>  (60) 

a 

Finally,  Clay  and  Tolstoy  show  that  the  scattering  signal 
can  be  expressed  as: 


4yza2>>l 


where 


<  S2  >  =  <pp*>  =  <PiPl*>  S, 

d  2  hf  . 

2  (2y2<C12>)  >>~ 


(61) 


A  =  ensonified  area 

<PiP?>!  the  expected  average  value  of  pt2,  where 
Pi  is  the  incoming  pressure  to  the 
luminated  area. 

Shf  =  scattering  function 

x2 

_ f2  (9)R2 _ 


e  (2yz<c 1 2>) 


2tt(cos  0i  +  cos  9Z)  <C  > 
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3 


1 


f (0)  =  f(0l0203) 


l+cos0  i  cos02-sinP  4  sin  9  2  cos 6  3 


COS0 1+COS0  2 

R  =  P ‘c ‘COSP  -  PC  COS01.  c  _  c 1 
p^^osG  +  pc  cos 0 1  sin0  sine1 

Jr 

Y  =  (cos  0!+  cos  0  2 ) 

tc  = 

COS  0 
J^ 

a  =  j  (sin  0!-  sin  02  cos  03) 


<^l2>  =  (3  +  5.12W)  x  10  W  —  wind  speed  in  m/s. 

Equation  (61)  is  valid  for: 

4y2a2>>l 

(2y2<;l2>)>> 

In  summary,  the  reflection  of  high  frequency  signals  yields 
scattered  radiation  which  is  incoherent.  Furthermore,  as 
pointed  out  by  Clay  and  Tolstoy  in  [Ref.  15]  although  the 
radiation  is  primarily  scattered  in  the  specular  direction, 
parts  are  scattered  in  all  directions.  As  seen  from 
Eq.  (61)  the  scattering  function  is  primarily  dependent 
on  the  mean  square  slope  of  the  surface  <<;12>  and  neither 
the  mean  square  wave  height  a2  nor  the  correlation  distance 
are  important.  Furthermore,  it  is  noticeable  that  in  the 
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high  frequency  limit  the  scattering  function  is  independent 
of  the  frequency  since 


e 


k2 

( 2 y  2 < ^  1  2  >) 


-*■1  when  r2-*-« 


C.  GEOMETRICAL  SHADOWING 


The  phenomena  of  shadowing  of  some  surface  areas  by 
others  has  to  be  considered  either  when  the  surface  irreg¬ 
ularities  are  large  with  respect  to  the  incident  wavelength 
or  when  the  grazing  angle  is  small. 

The  few  papers  devoted  to  this  subject  are  aimed  mainly 
towards  calculation  of  a  "Shadowing  function"  based  on 
the  statistics  of  the  surface. 

An  explicit  method,  geometrical  shadowing,  has  been 
introduced  by  P.  Bechmann  [Ref.  17]  where  the  shadowing 
function  S(9)  is  the  probability  that  the  point  ?  (Fig.  25) 
is  illuminated. 


S  ( 0 ) 


exp 


[-  j q  (x)  dx] 
o 


(62) 


where : 


q(x)  is  the  probability  that  c  is  shaded  by  <;  in 
the  interval  (x,x+dx)  given  that  it  is  not  shaded 
in  (0 ,x) . 

This  calculation  of  S(9)  only  considers  the  elevation  of 
the  surface  observation  point.  However,  the  slope  also 
plays  a  role  in  that  if  its  value  exceeds  cot8  the  point 
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will  certainly  be  shaded.  Thus,  R.  Wagner,  in  Ref.  18, 
incorporated  both  5  and  t, 1  using  P.  Bechmann's  method  and 
found  that  the  conditional  probability  that  a  point  on  the 
surface  is  illuminated,  given  that  it  has  height  5  and 

slope  c1/  can  be  expressed  as: 

00 

S(8,|?i/?|1)  =  exp [-  J  q(x)  dx]  u(cote  -  £ 1 )  (63) 

a 

where 

u:  is  the  unit  step  function. 

q(x):  is  the  conditional  probability  that  t,  is 

shadowed  in  the  interval  (x,x+dx)  given  that 
it  is  not  shadowed  in  (0,x) . 

The  function  q(x)  cannot  be  calculated  exactly.  Thus, 
R.  Wagner  made  the  approximation  that,  for  all  x,  the 
probability  that  ;  crosses  the  ray  in  dx  is  independent  of 
the  values  of  ?  and  s1  at  x=Q.  In  the  above,  no  mention 
has  been  made  of  the  direction  of  observation.  However, 
in  this  respect,  R.  Wagner  [Ref.  18]  points  out  that  in 
the  high  frequency  limit  only  those  portions  of  the 
surface  which  are  illuminated  simultaneously  by  rays  in 
the  direction  of  incidence  and  observation  can  contribute 
to  the  observed  scattered  power.  For  this  condition, 

R.  Wagner  introduced  both  directions  as  independent 
variables  in  the  so-called  bistatic  shadowing  probabilities 
where  he  defines: 
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1.  S(9i02lCi#^i')  to  be  the  conditional  probability 
that  the  surface  will  not  cross  the  incoming  ray 
(Ray  1)  or  the  outgoing  ray  (Ray  2)  anywhere, 
given  that  both  rays  pass  through  an  arbitrary 

point  or  the  surface  with  displacement  q  and  slope  ? 1 . 

2.  S(9i [e2, Si, Si')  to  be  the  conditional  probability 
that  the  surface  does  not  cross  Ray  1,  given  that 
it  does  not  cross  Ray  2  and  that  both  rays  pass 
through  the  point  ?  having  slope  c 1 . 

Thus,  the  conditional  shadowing  function  can  be  expressed 


as: 


S(0i,02Ui,Ci')  =  S(e1fe2*Ci»Ci*)  S(02|Ci#Ci’)  (64) 

The  shadowing  function  is  then  obtained  by  averaging 
over  all  possible  heights  and  slopes 


S  (6  i ,  e  2 ) 


oo 

jf  S(0,, 


e2U,?')  w(c,c')dc  dc' 


(65) 


Here  w(c,?')  is  the  bivariate  PDF  of  the  surface  height 
and  slope,  assumed  to  be  Gaussian 


««.?•>  -  h  ejtp{‘^r  '2Tf7r,  ,66) 


where  *o2  and  4>' '  are  the  values  at  t*0  of  the  correla¬ 
tion  function  and  its  second  derivatives,  respectively. 

For  the  region  0<9<j  where  the  probability  of  crossing 
one  ray  is  assumed  independent  of  that  of  crossing  the 
other,  R.  Wagner  found  that  the  bistatic  shadowing  function 
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feMbV'- 


as  we  from  Eq.  (60)  have  that  a2|t"(0)j 
=  <?12>  «  S2 

n i  =  cot  9;  i  =  1,2 
and  noting  that 

2  -t 2 

Error  function  erf(x)  *=  —  /  e  dt 

✓?  I 

2  00  ? 
f  —  t  2 

complementary  error  junction  ericixi  =  —  /  e  dt 

x 

The  shadowing  function  S(0i,02)  is,  in  short,  the  fraction 
of  the  surface  still  illuminated.  As  seen  from  Eq.  (60) , 
the  scattered  field,  in  the  high  frequency  case,  is  pro¬ 
portional  to  the  illuminated  area.  Hence,  the  shadowing 
effect  of  a  rough  surface  can  be  introduced  by  multiplying 
the  ensonified  area  A  by  the  shadowing  function  S(9l,82). 
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D.  ESTIMATING  THE  SURFACE  SCATTERING  EFFECT 


The  following  estimations  are  based  on  calculations 
in  the  specular  direction,  which  is,  as  pointed  out 
earlier,  expected  to  give  the  maximum  supplementary 
scattering  effect.  Hence,  in  the  specular  direction  where 
0!  =  3 ^  =  9  and  83  =  0°,  the  scattering  function  S^ 
reduces  to  the  following  expression: 


f2  (8)  R2 _  ( 

2  (cosG i+cos0 2 ) 2<s 1 2 > 


s  =  _  Rl 

h£  8tt<  C  1  2> 

as 

f ( 0 )  =  cos0 
Y  =  kcos0 
a  =  0 

x  =  —2—  =  0 
COS0 

As  pointed  out  in  the  previous  paragraph,  the  reflection 
of  very  high  frequency  signals  by  the  sea  surface  yields 
scattered  radiation  that  is  incoherent  under  the  assumption 
that 

4yza2  >>1 
(2y2<?1  2>)  >>^5- 

Before  we  launch  into  the  calculations,  we  will  verify 
these  criteria  for  the  frequency  range  of  interest:  60  and 
30  kHz. 
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Utilizing 


Y2  =  k2  cos20 


a  =  8. 1  x  10 


aW4 

46g2 


P  =  0.74 


<C12>  =  (3  +  5.12W)  x  10-3  g  =  9.81  m/s2 
and  assuming  a  windspeed  of  10  m/s  (SS3)  yields 
a2  =  0.30  m2  ;  a  =  0.55  m 


<C12>  =  5.42  x  10“2 

For  60  kHz(X  =  .025  m)  and  low  grazing  angles,  e.g.,  0  =  85° 
4y2o2  =  5.75  x  102  >>1 

(2y2  <£ 1 2>)  =  5.2  x  101  >>  - 

R2 

and  the  first  order  5.2>>—  ,  R^0.5m. 

R2 

Thus,  for  the  60  kHz  case,  the  criteria  are  fulfilled. 

For  30  kHz  (X=.05m)  and  9=85°. 

4a2 y 2  =  1.44  x  102>>1 

(2y2<£12>)  =  13>>  \ 

R 

again  the  first  order  R>1.0  m. 

Thus,  also  for  the  30  kHz  case  the  criteria  are  fulfilled. 

Similarly,  for  the  shadowing  function  in  the  specular 
direction 

O<9i*02»9<  rr/2  and  0  3  “  0 
we  obtain  the  following  simplified  expression 
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(69) 


S  (9)  =  S(01#02)  =  l1,-gx|J-4B)jerfv 


as 


Vi  =  V2  =  V  = 


(2a2  |  i)i  "  | 15  (2>2) 


>2^ 


n  1  —  n  2  =  n  =  cote 


B 1  ~  B  2  —  B 


[exp  (-v2  ? -/F  v  erfc  v] 


4  /rr  v 

In  summary,  for  specular  scattering  at  the  high  frequency 
limit,  the  expected  average  value  of  p2 ,  where  p  is  the 
pressure  field  at  the  receiver,  is  then  obtained  from  the 
following  simplified  expressions: 


where 


<s2>  *  <pp*>  s 

Rz 


<S*>  =  <pp.>  “1515! 


hf 


3ttR2  <C12> 


S ( 0 )  _  1-exp (-4B)erfv 


4B 


(70) 


v  = 


( 21 2 ) 


B 


COt0 

[exp  (-v2 )  -/rfv  erfc  vl 
4  (  tt  )  ^  v 
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<£12>  =  £2  =  (3  +  5.12W)  x  10  ^  ;  W  =  windspeed  in 
m/s. 


R  -  p1c1cos9  -pc  cos6l 


P^^osO  +  P  c  cosQ1  sin®  sine1 


From  Ref.  19  we  use  the  following  air/sea  water  inter¬ 
face  data: 

1.  For  air  c'  =  343  m/s  and  p^1  =  415  Ray  Is. 

2.  For  sea  water  c  =  1500  m/s  and  pc  =  1.54  x  106  Rayls 
The  estimation  of  the  illuminated  area  A  for  specular 

scattering  where  0i  =  02  =  0  and  03  =  0,  is  based  on  the 
geometry  illustrated  in  Fig.  24b.  The  illuminated  area 
is  given  by 

A  =  -nab  (71) 

Assuming  the  following  data  to  be  known 
hs  =  depth  of  the  source 

h  =  depth  of  the  receiver 
r 

Atp  =  the  half  beam  width  of  the  directional  receiver 
Both  a  and  b  of  Eq.  (71)  can  be  calculated  in  terms  of  the 
detection  range  R_  as  follows: 


*2 


*D  hr 
h  +h 


(72) 


Xi 


«D  * 


X  2 


V1 


h 

■e“tk: 


tan-1 ( 


*0 


V*r 
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Ri  = 


sin9 


R  2  = 


sxn0 


To  a  good  approximation  when  9  =  ir/2 ,  we  have 


a  -  x.  -  V1  -  ETT1 

s  r 


b  =  R  sin  (A<j>)  = 


»D  hr 


(hg+TTj  sm0 


sin  ( A4> ) 


The  expression  for  <s2>  is  then  introduced  as  a  supple¬ 
ment  to  the  direct  path  to  the  receiver  in  the  following 


way: 

By  utilizing  the  relationship 

I  =  E  (74, 

for  the  intensity,  the  scattering  intensity  at  the  receiver 
and  the  intensity  at  the  ensonified  area  are,  respectively, 

i  =  E  =  <s-.>- 

s  pc  c 


P2  <P1P1*> 

;  .  a  — -  =  ■  —  ■  — <  ■ 


1  Pc  PC 

Equation  (70)  can  thus  be  written  as 

I  *  t  AS  (9)  R2 
3  1  8ttR2  2<c  1 2> 


V1* 


AS (9)  R2 
8trR2  2<<; 1  2> 


•ft- <• 


10 log  Ig/I, 


lOlog 


I  ,1  - 

s/  ref  _ 


ref 


K 

o 


IL  =  IL (R , )  =  K 
s  1  o 


(75) 


We  then  have  to  determine  IL(Rj) 


TL  =  SPL(l)  -  SPL(Rj)  =  201ogR1 
SPLfRj)  =  IL (R  t )  =  lOlogl t 
SPL(l)  =  SL 

yielding 

IL(R1)  *  SL  -  201ogRt  =  Kx 

Then 


ILsuppl  *  Ko  + 

is  the  supplement  to  the  direct  path,  and 
ILdirect  *  SL  -  J01o«bd  *  K> 

W'raJ  ‘  anti  l°?  ’  K- 

ILs 

Isuppl/1re£  *  anti  1o?  UT  *  K’ 

The  total  intensity  of  the  receiver  is  thus 


ITot/Iref 


•^direct  +  Isuppl 
^■ref 


K„  +  Ks 


(76) 


(77) 

(78) 


ILTQt  =  101ogK6  =*  K7 


(79) 
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Thus,  the  effect  of  the  randomly  rough  surface  compared 
to  the  idealized  free-field  condition  can  be  expressed  as 


AIL  =  IL_  .  -  IL, .  ,  =  Ks 

Tot  direct  8 


(80) 


A  calculator  program  on  a  Texas  Instrument  59  (later 
called  TI  59)  was  developed  to  perform  these  calculations. 
A  block  diagram  of  the  program  is  outlined  in  Fig.  (26) 
and  the  programs  steps  together  with  a  detailed  description 
is  given  in  Appendix  B.  The  calculations  are  based  on  the 
following  fixed  data 
h  =  2  m 


s 

hr  =  6  m 
AcJ>  -  10° 

W  =  10  m/s  (SS3) 

Then,  varying  the  detection  range  from  R  =  2000  m  to  100  m 
gives  the  difference  between  and  ILcjirect  plotted  in 

Fig.  (27). 


E.  CONCLUSION  AND  DISCUSSION 

As  seen  from  Fig.  (27)  the  supplementary  effect  of 
the  scattering  from  a  rough  surface  in  the  high  frequency 
case  is  negligible  compared  to  the  direct  path. 

In  saying  so,  it  also  should  be  pointed  out  that  the 
Helmholtz-Kirchhof f ' s  approach  may  be  limited  as  it  does 
not  take  into  account  the  diffraction  effects  from  crests 
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and  throughs  of  the  ensonified  area,  an  effect  which 
becomes  increasingly  important  at  low  grazing  angles 
high  frequency  and  when  the  rough  surface  is  a  super 
position  of  swell  and  capillary  waves. 
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VII.  THE  EFFECT  OF  SCATTERING  AND  ABSORPTION  FROM  THE  SUB- 
SURFACE  OCEAN  LAYER  ‘ 


A.  GENERAL  SCATTERING  THEORY 


Generally  when  a  region  (volume)  scatters  sound,  some 
of  the  energy  carried  by  the  incident  wave  is  dispersed. 

The  energy  lost  by  the  incident  wave  may  be  absorbed 
by  the  scatterers  or  it  may  be  simply  deflected  from  its 
original  course.  The  amount  of  energy  lost  per  second  by 
the  incident  wave  divided  by  the  incident  wave's  intensity 
is  called  the  total  cross  section  a g  of  the  region  and  is 
the  sum  of  the  absorption  and  scattering  cross  sections 


where 


a 

e 


n  +n 

s 

I 


p 


a 


ng  =  scattered  power 


n  =  absorbed  power 

cl 


The  existence  of  gas  bubbles  in  the  subsurface  ocean 
layer  modifies  the  forward  scattering  in  the  following 
two  major  ways: 

1.  The  bubbles  can  resonate.  When  the  bubbles  are 

excited  at  a  frequency  near  its  natural  frequency, 
it  very  efficiently  absorbs  and  scatters  the  inci¬ 
dent  wave.  At'  resonance,  the  scattering  and 
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I 


I 

i 


I 


I 


r 


i 


absorption  cross  section  of  a  typical  bubble  at  sea 

3 

is  of  the  order  10  times  its  geometrical  cross 
section. 

2.  The  bubbles  change  the  effective  compressibility  of 
the  water  and  cause  the  speed  of  sound  to  be  a 
function  of  frequency,  i.e.,  the  medium  is  dispersive. 

We  will  investigate  and  discuss  these  effects  by  separately 
estimating: 

1.  The  attenuation  due  to  the  bubbles  and 

2.  the  refraction  by  bubbles. 

B.  ABSORPTION  MODEL 

The  choice  of  model  for  the  subsurface  ocean  layer  depends 
on  whether  the  medium  has  a  teneous  or  a  dense  distribution 
of  scatterers.  When  the  bubble  density  is  teneous,  both 
"single  scattering"  and  "first  order  multiple  scattering” 
approximation  solutions  are  applicable. 

On  the  other  extreme,  when  the  bubble  density  is  high, 
the  so-called  "diffusion"  approximation  can  be  used. 

Between  these  two  extremes,  multiple  scattering  effects  are 
important . 

The  multiple  scattering  theory,  which  in  the  limit  also 
contains  the  first  order  approximation,  will  be  used  to 
estimate  effects  of  attenuation  due  to  bubbles,  on  the 
propagation  of  propeller  noise  from  the  target  to  the 
torpedo . 


1 
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The  geometry  of  the  propagation  model  is  illustrated 
in  Fig.  (28)  where  it  is  assumed  that  a  plane  wave  is 
incident  on  a  semi-infinite  (disregarding  the  sea  surface) 
slab  of  thickness  x  containing  a  number  of  randomly  dis¬ 
tributed  bubbles.  The  plane  wave  approximation  is  valid 
if  the  incident  sound  has  a  wavelength  X  much  greater 
than  the  bubble-radius  a 
ka<<l 

where 

k=(jj/c=2T/X 

The  receiver  is  located  outside  the  slab  and  the  beam 
pattern  of  the  receiver  is  represented  by  the  solid  angle 

V 

We  are  interested  in  the  estimating  of  the  total  power 
received,  taking  into  account  the  multiple  scattering 
process  in  the  inhomogeneous  slab  as  well  as  the  beam 
pattern  of  the  receiver. 

The  mathematical  formulation  of  this  problem  is  based 
on  Twersky's  theory  of  multiple  scattering.  Since  the 
theory  is  presented  in  Ref.  20  only,  the  basic  formulation, 
major  assumptions,  and  the  end  results  will  be  presented 
here. 

The  total  intensity  is  the  average  of  the  square  of 
the  magnitude  of  the  total  field: 
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<l^alz>  =  <!<<i>a>  +  i^fai2>  =  |«j>a>|2  +  <Ufa|2>  (8i) 

where 

a  a  N 

1)'  -  'K  +  Z  u  a 

1  s-1  3 

the  scalar  field  at  the  receiver  location  r  , 

cl 

see  Fig.  (29)  ,  is  the  sum  of  the  incident  wave  4>^ 
and  the  contribution  from  all  N  scatterers. 

|  < i^a>  |  2  is  the  coherent  intensity  based  on  the 
average  field  <  tfia> .  <|i|;^a|2>  is  the  incoherent 
intensity  based  on  the  fluctuating  field  ^^a. 

In  Twersky's  theory,  the  multiple  scattering  process  is 
described  by  the  following  set  of  integral  equations  which 
Eq.  (81)  must  satisfy: 

<|'(/a!2>  -  |<<(/a>|2  +  |  v  a|  2  <|^s|  2>p(r  )dr  (82) 

S  5  S 


where 


*.*  •  usa  + 


Uta  vsfc  p(rfc)dt 


is  an  operator  representing  all  the  scattering  processes 
from  s  to  a.  (See  Fig.  (28).)  It  should  be  noted  that 
Twersky's  theory  includes  all  the  multiply  scattered  waves 
that  involve  chains  of  successive  scattering  going  through 
different  scatterers.  (See  Fig.  (29a).)  However,  the  theory 
neglects  the  terms  which  include  an  individual 
scatterer  more  than  once,  as  illustrated  by  Fig.  (29b) . 
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Thus  Twersky's  theory  is  expected  to  give  good  results 
when  back  scattering  is  insignificant  compared  to  the 
scattering  in  other  directions. 

As  typical  for  most  integral  equations,  Twersky  being 
no  exception,  detailed  complete  solutions  are  difficult 
to  obtain.  However,  Twersky  gave  an  approximate  solution 
to  Eq.  (81)  and  Eq.  (82) ,  which  according  to  Ref.  20  have 
been  found  to  agree  reasonably  well  with  experimental  data. 
This  solution  is  based  on  the  following  main  assumptions: 

-Backscattering  is  assumed  to  be  small  compared  to 
scattering  in  other  directions. 

-Scattering  is  mostly  concentrated  in  the  forward 
direction.  This  is  reasonable  based  on  the  assumed 
random  distribution  of  the  bubbles;  i.e.,  no  rein¬ 
forcement  of  the  radiation  pattern  occurs  except  in 
the  direction  of  the  incident  wave. 

-The  angle  9  is  small,  i.e.,  0  =  0°. 

u  S  a  5 

This  leads  to  the  following  expression  for  the  total 
intensity  at  the  receiver: 

<l^a|2>  *  T  =  exp(-pa  x)  exp (-pa  x)  (83 

cl  S 

+q[l-exp(-pagx) ] 

where 

a  58  absorption  cross  section 

3l 

a  *  scattering  cross  section 
s 
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..4|f|Ma» 


4  TT 


f\f\2dn 

TT  '  1  C 


where 

f  =  the  amplitude  function 

q  =  the  fraction  of  total  scattered  power  collected  by 
the  receiver  as  illustrated  in  Fig.  30  and 

p  =  is  the  bubble  density,  i.e.,  the  number  of 
scatterers  per  unit  volume. 

For  small  values  of  px  we  see  from  Eq.  (82)  that  the 
coherent  part  dominates: 

InT  =  -(a  +a)px  (84) 

a  S 

In  this  limit  the  multiple  scattering  result  is  equiva¬ 
lent  to  that  obtained  from  single  scattering  considerations. 

For  large  values  of  px  which  corresponds  to  very  dense 
or  very  wide  slab  of  scatterers,  the  incoherent  intensity 
dominates 

InT  ~  lnq  -  a  px  (85) 

d 

In  this  latter  case,  it  is  notable  that  when  Qr=2ir , 
the  receiver  collects  almost  all  the  scattered  power 

InT  :  -  o  px  (86) 

d 

The  first  case,  representing  the  situation  for  teneous 
density  of  scatterers  and/or  narrow  beam  pattern  of  the 
receiver,  gives  a  good  approximation  to  the  situation  of 
interest  in  the  thesis.  It  also  represents  the  case  for 


74 


which  no  scattered  power  is  received.  This  will  be 
approximately  true  for  a  narrow-beam  width  receiver.  We 
will,  therefore,  investigate  the  coherent  intensity  first. 

Also,  the  incoherent  case,  as  represented  by  Eq.  (86) 
will  be  investigated,  where  only  losses  due  to  absorption 
are  incorporated. 


C.  THE  COHERENT  INTENSITY  CASE 

For  the  coherent  case  the  intensity  level  after  the 
incident  wave  has  traversed  a  distance  x  is 


where 


Ix  *  Ip  exp[-(aa+cg) px] 


Ip  =  incident  plane  wave  intensity. 


The  change  in  intensity  over  the  distance  x  is 


AIL  =  lOlog 


!x  _  Ip  exp  [-(aa+Ps)px] 


(87) 


(88) 


=  lOlogexp  [-(a  +a _)  px] 

a  S 

The  excess  attenuation  per  unit  distance  due  to  bubbles 
is  thus 

- (o  +a)  px  log  e 


AIL  =  _  r 
x  **  *•' 


(°  +a  )  p log  e 

d  S 


a  s 


-] 


a  *  4.34  a  p,  in  dB/m. 


(89) 


However,  this  only  takes  into  account  bubbles  of  one 
size.  In  a  bubbly  medium  there  is  a  spectrum  of  radii. 


The  probability  density  function  for  finding  a  bubble  size 
between  radii  a  and  a+da  is 

W  (a)  =  2^-  (90) 

where 

00 

/  W  (a)  da  =  1 

o 

n(a)da  is  the  number  of  bubbles  per  unit  volume 
having  radii  between  a  and  a+da.  It  is  common  to 
use  da=l  urn. 

CD 

p  =  n(a)da  is  the  total  number  of  bubbles  per 

o 

unit  volume. 

As  the  extinction  cross  section  also  is  a  function  of 
the  radius  (See  Eq.  (93).),  the  absorption  due  to 
bubbles  is  obtained  by  integrating  Eq.  (89)  over  all 
possible  radii 

00 

a  =  4.34  /  ae(a)n(a)da  (91) 

o 

To  calculate  the  absorption  coefficient,  the  extinction 
cross  section  ag  must  be  derived  from  the  general  bubble 
dynamic  relationship.  This  is  done  in  detail  by  C.  Clay 
and  H.  Medwin  [Ref.  21]  from  which  the  following  results 
are  taken. 
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The  scattering  cross  section 


where 


°S  -  r 


[(fj/f)  2-l]  2  +  62 


2ybgp 


1  M  A 

fr  =  resonance  frequency  =  ( - )  3 

f  -  f  =  operating  frequency 
6  *  damping  constant  =  5  +  <5fc  +  6v 

=  ka  +  (|)  (^)2  +  r^rr 


a  *  bubble  radii 


PAu>a< 


d  _  ^  (y_d  [H(sin  h  X  +  sin  X)  -2  (cos  h  £  -  cos  £) _ ^ 

^  X2 (cos  h  X  -  cos  X)  +  3(y-l) (sin  h  X  -  sin  X 

3w  C  , 

x  =  a(_P3^a)'i 


kg  =  thermal  conductivity  of  gas 

pg  =  density  of  gas  =  PgAU  +  ^p  a)  ]  (1+0. lz) 

A 

p  -  *  density  of  gas  at  sea  level 
gA 


t  »  surface  tension 
>  »  1.013  x  106 (1+0. lz) 


z  *  bubble  depth  in  m. 

Cpg  *  specific  heat  of  constant  pressure  of  gas 


U  =  shear  viscosity  of  water 
Y  =  7/5,  for  diatomic  gas 


0^  =  density  for  sea  water 


b  =  [l+(|)*r‘  [1  +  sin  h JC  -  sinjj 

'J  X  cos  h  X  -  cos  X 


(93) 


(94) 


A  detailed  computer  program,  as  outlined  in  Appendix  B, 
was  developed  for  the  TI  59  to  handle  the  derivation  of 
u  .  o  ,  and  a  based  on  an  assumed  receiver  depth  of  z  =  6  m. 

For  the  60  kHz  case,  both  the  extinction  cross  section 
o  and  the  absorption  cross  section  o  are  given  in  Fig.  31 
as  a  function  of  bubble  radius  a. 


Similarly,  Fig.  32  gives  o  and  a  for  the  30  kHz  case. 

6  cl 

Superimposed  on  these  figures  are  the  curves  for  n(a)da  as 
calculated  from  the  following:  Figures  5  -  7  of  Ref.  2 
give  the  resonant  bubble  densities  in  a  1  urn  band  as  a 
function  of  depth  and  with  the  wind  speed  as  parameter  for 
the  three  discrete  frequencies  12  kHz,  38  kHz  and  120  kHz. 
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Based  on  these  data.  Fig.  33  shows  the  interpolated  bubble 
density  as  a  function  of  resonance  frequency  for  sea  state 
2 ,  3 ,  and  6 . 

Furthermore,  both  A.  Lovik  [Ref.  2]  and  H.  Medwin  [Ref.  3] 

found  that  the  bubble  density  function  n(a)  decreases  with 

increasing  bubble  radii  as 
-x 

naa 

where  H.  Medwin  [Ref.  3]  found  the  power  law  of: 
x  =  4  for  a<50-80  ym. 
x  =  2  for  a>50-80  ym. 

and  A.  Lovik  [Ref.  2]  found  the  power  law  of: 
x  =  4.2  between  38  kHz  and  120  kHz 
x  *  2.6  between  12  kHz  and  38  kHz 
which,  averaged  over  the  depth  interval,  corresponds  to  the 
bubble  radii  of  380  ym(12  kHz),  120  ym  (38  kHz)  and  49  ym 
(120  kHz) . 

As  suggested  by  A.  Lovik  [Ref.  2],  the  discrepancy 
between  the  two  observations  is  not  great  and  may  be  due 
to  the  few  measuring  frequencies  used  in  the  work  of 
A.  Lovik  . 

In  summary,  Fig.  33,  from  which  we  obtain  the  appro¬ 
priate  resonant  bubble  density  in  a  1  ym  band  n(aR) 
together  with  the  power  law  naa~x,  comprise  the  full 
knowledge  of  n ( a ) da . 
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Performing  a  multiplication  of  aQ  and  n(a)da,  we 
obtain  Fig.  34  and  Fig.  35  for  the  60  kHz  and  the  30  kHz 
case ,  respectively . 

Finally,  the  integral 

CO 

/  oe(a)n(a)da 

o 

was  evaluated  using  a  numerical  integration  based  on 
Simpson's  discrete  approximation  programmed  for  the  TI  59 
and  documented  in  Appendix  C. 

Based  on  the  above,  the  following  absorption  coeffi¬ 
cient  for  the  coherent  case  is  obtained  for  f  =  60  kHz 

CO 

a  =  4.34  /  a  (a) n (a) da  =  4.34  (2.016  x  10_1) 
o  e 

■  8.75  x  10_1  dB/m. 

a  *  0.88  dB/m  (95) 

with 

n (aR)  *  1000 
n(a)  a  a-4 
z  =  6  m. 

Thus,  it  is  seen  that  the  attenuation  due  to  bubbles 
is  considerably  greater  than  the  "normal  attenuation" 
due  to  chemical  and  viscous  relaxation  processes  in  sea 
water,  which  for  the  60  kHz  case  is 
a  ■  0.02  dB/m. 
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Thus,  the  total  absorption  coefficient  for  the  60  kHz 
case  is 

a  =  0.90  dB/m.  (96) 

For  f  =  30  kHz,  the  absorption  coefficient  due  to  bubbles 
in  the  coherent  case  is 

oo 

a  =  4.34  /  ag (a) n (a) da  =  4.34  (3.8375  x  10-2) 

o 

=  1.665  x  10_1  dB/m. 

a  ~  0.17  dB/m  (97) 

with 

n(aR)  =  20 
n  (a)  a  a  2  ’  ® 
z  =  6  m. 

-2 

For  this  case,  the  power  law  dependence  of  n(a)  oc  a 
gives 

a  *  1.655  x  10_1 

Thus,  the  difference  in  power  law  dependence  makes  no 
significant  difference  in  the  absorption  coefficient. 

The  absorption  coefficient  due  to  chemical  and  viscous 
relaxation  processes  is  at  30  kHz 
a  =  0.012  dB/m. 

The  total  absorption  coefficient  in  the  coherent 
case  is 

a  =  0.18  dB/m.  (98) 
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D.  THE  INCOHERENT  INTENSITY  CASE 

For  the  incoherent  case,  where  only  losses  due  to  ab¬ 
sorption  are  included,  the  product  a  and  n(a)da  for  the 

cl 

60  kHz  and  30  kHz  cases  are  given  in  Fig.  36  and  Fig.  37, 
respectively.  Performing  a  numerical  integration  based  on 
the  Simpson's  discrete  approximation  leads  to  the  follow¬ 
ing  results : 

For  f  =  60  kHz,  the  absorption  coefficient  due  to 
bubbles  is 

GO 

a  =  4.34  /  a  (a)n(a)da  =  4.34  (1.66  x  10-1) 

0 

a  s  0.72  dB/m.  (99) 

with 

n(aR)  =  1000 

,  x  -4 
n(a)  a  a 

z  =  6  m. 

Adding  the  "normal  attenuation"  in  sea  water  for 
f  =  kHz,  yields  a  total  absorption  coefficient  of 

a  *  0.74  dB/m  (100) 

For  f  =  30  kHz,  the  absorption  coefficient  due  to 
bubbles  is 

CO 

a  =  4.34/a  (a)  n  (a)  da  =  4.34  (2.975  x  lO-2)  (101) 

0  a 

a  =  0.13  dB/m 


82 


with 


n(aR)  =  20 

n(a)  a  a  ^ 
z  =  6  m. 

Adding  the  "normal  attenuation”  in  sea  water  for  f  =  30  kHz 
yields  a  total  absorption  coefficient  in  the  incoherent  case 
of 

a  s  0.14  dB/m  (102) 

E.  SUMMARY  AND  DISCUSSION  OF  THE  BUBBLE  ATTENUATION 

As  seen  from  the  above,  the  attenuation  due  to  bubbles  in 
the  subsurface  ocean  layer  is  important  for  high  frequency 
and  high  Sea  State. 

Accounting  for  the  effect  of  bubbles  at  Sea  State  3  in 
summary  we  found  the  following  absorption  coefficients  in 
dB/m: 

f  =  60  kHz  f  *  30  kHz 

coherent  incoherent  coherent  incoherent 

a  =  0.88  a  =  0.72  a  =  0.17  a  =  0.13 

The  main  assumptions  were: 

-the  back  scattering  is  small  compared  to  that  in 
other  directions. 

-the  scattering  is  mostly  concentrated  in  the  » 

* 

forward  direction. 


i 
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VC*. 


-the  angle  0  between  the  scatterer  and  the 

3.S 

receiver  is  small,  i.e.,  0  =0°. 

as 

-the  depth  of  the  receiver  is  2  =  6  m. 

-the  sea  state  is  3. 

The  range  dependent  portion  of  the  passive  sonar 
Eq.  (1)  TL=-201ogR-aR  for  both  the  60  kHz  and  30  kHz 
scattering  results  can  now  be  compared  with  the  reference 
data  obtained  in  Section  V. 

For  f  =  60  kHz,  Fig.  22  gives  TL=-201ogR-aR  as  a 
function  of  R  with 

1.  a  =  0.021  dB/m,  the  "normal  attenuation”  due  to 
chemical  and  viscous  relaxation  processes. 

2.  a  =  0.9  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  coherent  case. 

3.  a  =  0.74  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  incoherent  case. 

For  the  same  example  as  in  the  reference  model,  Fig.  22 
yields  the  detection  ranges  for  60  kHz. 

R  =  1200  m  for  a  =  0.021  dB/m  (103) 

R  =  60  m  for  a  *  0.90  dB/m 
R  =  70  m  for  a  =  0.74  dB/m 

Not  surprisingly,  this  result  seems  to  exclude  the  possi¬ 
bility  or  having  both  a  searching  and  attack  depth  near 
the  surface,  i.e.,  z  =  6  m,  for  a  torpedo  system  operating 
in  a  high  frequency  region,  f  *  60  kHz. 
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Similarly,  for  f  =  30  kHz,  Fig.  23  gives  TL=-201ogR-aR 
as  a  function  of  R  with 

1.  a  =  0.012  dB/m,  the  "normal  attenuation"  in  sea 
water. 

2.  a  =  0.18  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  coherent  case. 

3.  a  =  0.14  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  incoherent  case. 

For  the  same  detection  example  as  in  the  reference  model. 
Fig.  23  yields  the  detection  ranges  for  30  kHz: 

R  =  2400  m  for  a  =  0.012  dB/m  (104) 

R  =  250  m  for  a  =  0.18  dB/m 

R  =  310  m  for  a  =  0.14  dB/m. 

Again,  the  bubbles  give  a  major  decrease  in  the  detection 
range.  A  detection  range  of  R  =  250  m  seems  marginally 
acceptable  as  the  turn  rate  requirement  for  a  pursuit 
homing  trajectory  may  become  excessive. 

The  above  results  are  summarized  in  Tables  III  and  IV 
for  the  60  kHz  and  30  kHz,  respectively. 

F.  THE  REFRACTION  BY  BUBBLES 

The  presence  of  bubbles  in  the  sea  water  affects  the 
speed  of  sound  (phase  speed)  primarily  because  of  the  change 
in  compressibility.  The  derivation  of  this  dispersive  rela¬ 
tionship  on  the  sound  speed  has  been  done  by  H.  Medwin 
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fc.  id  If  1 1  ■« 


!“Ref .  22]  .  H.  Medwin  showed  that  the  bubbles  with  resonant 
frequencies  greater  than  the  incident  frequency  decrease 
the  sound  speed,  while  bubbles  with  resonant  frequency  lower 
than  the  incident  wave  increase  the  sound  speed. 

Furthermore,  H.  Medwin  [Ref.  22]  predicts  the  sound  speed 
gradient  due  to  bubbles  as  a  function  of  depth  for  a  frequency 
range  and  wind  speed  compatible  with  our  domain  of  interest. 

He  found  the  gradients 

g  =  3c/3z  =  0.26  s  ^  at  z  =  0  m. 
g  =  3c/3z  =  0.016  s  ^  at  z  =  10  m. 
g  =  3c/3z  =  0.005  s”1  at  z  =  20  m. 

For  comparison,  the  sound  speed  gradient  due  to  pressure 
in  an  isothermal  water  is 
g  =  0.017  s"1 

This  shows  that  the  rays  in  the  top  10  m  are  strongly 
influenced  by  bubbles.  However,  with  respect  to  our  surface 
scattering  model  (Section  VI)  where  both  the  source  and  the 
receiver  are  situated  very  close  to  the  surface  and  where 
relatively  short  propagation  distances  are  encountered, 
this  refraction  phenomenon  is  assumed  to  have  negligible 


effect. 


VIII.  THE  TURN  RATE  LIMITATION 

As  pointed  out  in  the  previous  section,  the  presence 
of  bubbles  near  the  surface  may  significantly  reduce  the 
range  at  which  the  target  can  be  detected. 

In  this  section,  the  turn  rate  necessary  during 
pursuit  homing  at  the  previous  estimated  detection  ranges 
will  be  studied.  A  computational  procedure  will  be  used 
to  determine  the  range  of  angles  on  the  bow  (AOB)  of  the 
target  at  the  beginning  of  homing  which  lead  to  miss  on 
the  initial  attack,  A  trajectory  where  the  torpedo  velocity 
vector  always  is  directed  towards  the  instantaneous  target 
position  is  called  a  pursuit  homing  trajectory.  The 
derivation  of  the  pursuit  homing  trajectory  follows 
P.  van  Nostrand  [Ref.  23]  and  is  based  on  the  geometry  of 
Fig.  38,  where 

r  =  approach  angle,  i.e.,  angle  between  ship 
velocity  vector  and  the  line  of  sight. 

AOB  =  180  -  $ 

The  equation  of  motion  is  obtained  by  taking  the 
component  along  r  and  the  normal  to  r,  yielding 


r  =  V„  cos  <1>  -  Vm 

S  T 

(105) 

r<J  *  -V„  sin  $ 
s 

(106) 
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where 


r  =  range  rate 


=  turn  rate 


Vs  =  target  speed 


VT  =  torpedo  speed. 


Dividing  Eq.  (105)  by  Eq.  (106)  yields: 


I  -  (v;  ns? '  oot  01  *  ’  ri4° 


and  defining 


yields 


V_  p 


7  =  {ST7T*  "  cot  ♦ 
r  sm<j> 


/  •  *iP“l 
r-  -  T. (s3-n<j>) 

r  ~  *(l+cos<!>)  f 


(107) 


(108) 


Then,  introducing  the  initial  conditions:  rQ,  <pQ  where 

r  :  initial  detection  range 
o 


:  initial  approaching  angle 


yields 


K  =  r 


(1+cos  4>  ) 
°(sin  <t|0)P  ^ 


(109) 


For  the  geometry  of  Fig.  38,  the  turning  rate  is  given  by 
Eq.  (105). 
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(110) 


V 

I  =  sin  p 

Substituting  Eq.  (108)  into  Eq.  (110)  yields: 

P 


♦  =  - 


Vg(l+COS  0 ) 


(111) 


K (sin  4>) 

From  Eq.  (108)  we  see  that  r-*-0  as  <p-*0,  i.e.,  <j»  tends  to 

zero  as  the  torpedo  approaches  the  target  ship.  It  is 
further  of  interest  to  determine  the  limiting  value  of  the 
turning  rate  as  the  torpedo  approaches  the  target  for 
various  values  of  the  parameter  p.  This  is  done  by  taking 
the  derivative  of  <j>  with  respect  to  p.  Thus,  from  Eq.  (109) 
we  get 


-  -  -jf (sin  1-P  [l+cos<j>]P  [2  cos  p  -  p]  (112) 

For  l<p<2,  we  see  that  Eq.  (112)  is  zero  at 

cos  5>  =  P/2  (113) 

p  *  cos'1  (P/2) 
with  am  associated  maxima 


V  (l+P/2)P[l-(P/2) Z]1'P/2 
I ♦max I  "  K  (114) 

Furthermore,  the  turn  rate  for  p<2  at  impact  is  zero  as 

the  limit  of  Eq.  (109)  yields: 

lim  lim(sin  <J>)  2_P  [1+cos  <J>]P  -*-0  (115) 

4h-0  <t>-^0 
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For  p=2  we  see  from  Dq.  (Ill)  that 

V  ?  4V 

lim  0=  -  lim(l+  cos  <P)*  - - -  (116) 

0-0  “*"0 

For  p>2  we  see  from  Eq.  (Ill)  that  there  is  no  zero 
value  of  d0/d0  between  0  -  0°  and  0  =  180°,  since  all  terms 
of  Eq.  (112)  are  nonzero  terms  for  0°<p<180°.  Furthermore, 
the  turn  rate  for  p>2  at  impact  is 


I  lim  i  |  -  | -  ^  lim  (117) 

*.0  <sin 

as  (sin  0)-°°  since  (p-2)>0. 

Furthermore,  as  0-180°,  we  get  from  Eq.  (112) 

•  V 

lim  =  -^(2+p)  lim  (sin  0)1-P (1+cos  <j>)P 
0-180 

Evaluating 

lim  fcin  0)1_P  (1+cos  0)P  =  lim  (sin  p)  (jjS25-i) p 

sin  cp 

0—180  0-180 
by  applying  1'  Hopitale  rule  to  the  term 


1+cos  0  _  i _-sin  0 
lim  — -rr — z —  “  1  iin  ^ 

sin  0  cos  <p 

0-180  4>-180 


0 


1+cos  4)  P 

Thus,  the  product  sin  0 (  )  approaches  zero  for 

any  p>l,  since  both  terms  in  this  product  approach  zero 
as  0-180. 
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Figure  39  shows  a  plot  of  computed  values  of  j^K/V  j 

s 

for  different  values  of  p>l.  From  this  we  can  draw  the 
following  summarizing  conclusions. 

For  l<p<2  the  turn  rate  has 

-1 

-a  maximum  value  at  <t>=cos  (p/2) 

-the  zero  value  at  <j>=0°  and  0=180°. 

For  p=2  the  turn  rate 
-is  zero  at  <p=l 80°. 

-monotonically  increasing  with  decreasing  <p 
approaching  the  value  (4Vs/K)  when  cj>=0° . 

For  p  2  the  turn  rate 
-is  zero  at  $=180°. 

-is  monotonically  increasing  with  decreasing  $ 
approaching  <*>  as  0=0°. 

As  seen  from  the  above  for  p>2,  the  turn  rate  increases 
monotonically  with  decreasing  approaching  an  infinite 
turn  rate  to  hit  a  point  target.  To  avoid  this  singularly 
we  must  make  some  provision.  If  the  torpedo's  maximum 
turn  rate  is  exceeded  only  at  some  very  small  range,  a  hit 
is  likely. 

Figure  40  illustrates  the  hit  criterion  used.  Assume 
a  rudder  of  length  LR  is  situated  directly  behind  the 
propeller  which  is  idealizer  as  an  acoustic  point  source. 

If  the  torpedo  becomes  turn  rate  limited  at  some  range 
r=r',  it  at  best  can  proceed  along  a  circular  path  which 
lags  the  desired  pursuit  trajectory,  or  at  worst  it  can 
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loose  acoustic  contact  and  go  into  "hold-in,"  maintaining 
a  constant  heading  at  the  angle  4> '  .  In  this  last  case,  it 
will  cross  the  line  of  target  ship's  course  in  the  time: 

t  =  r'/V  (118) 

If  this  advance  is  less  than  the  length  of  the  rudder  LR, 
the  torpedo  will  impact  the  rudder.  Hence,  the  limiting 
condition  is 


r'  =  pLR  (119) 

For  this  analysis,  we  have  arbitrarily  chosen 

Lr  =  3  m  (120) 

Thus,  the  torpedo's  turn  rate  has  not  been  exceeded 
when  the  range  to  the  target  is 

r'  =  pLR  =  3  p  (121) 

and  a  hit  is  assumed. 

Now  we  can  analyze  the  cases  p>2  and  l<p<2  on  the  same 
footing.  If  we  match  the  torpedo's  maximum  turn  rate  to 
a  particular  range  r',  it  is  certain  that  its  turn  rate 
has  not  been  exceeded  earlier  in  the  run. 

This  turn  rate-range  matching  is  done  by  substituting 
the  "hit-range"  defined  as  r'  =  3  p  into  Eq.  (106)  and 
solving  it  for  sin  <p,  yielding 
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max 


Vs  .  4 

3p  * 


sin  0  =  - 


$  (3p) 

~r.ax  ^ 


Since  j>  is  always  negative,  sin  $  is  always  positive  and 
equal  to 

3p  f 


sin 


1  max 


V 


(122) 


Two  values  of  i  satisfies  this  equation,  and  they  are  denoted 

.  -1  3P|' 


<j>A  =  sin 


V 


for  0  ^  90° 


(123) 


*B  =  180  *  ^ 

A  unique  value  of  K  may  now  be  found  using  Eq.  (Ill)  for 
each  of  the  angles  $A  and  1>B. 


K  = 


Vs(l+cos  $) 


iL.v(sin  <P) 
‘  max 


P-2 


These  values  of  K  are  designated  as  KA  and  Kg,  and  for 
each  there  is  a  corresponding  value  of  <PQ  from  Eq.  (109) 
with  $  given  by  the  initial  detection  range 


(Va  and  (Vb 

In  the  case  of  l<p<2,  the  turn  rate  does  not  necessarily 
increase  monotonically  during  the  pursuit  homing  trajectory, 
and  we  must  check  that  either  $A  or  $B,  respectively,  is  not 


beyond  the  angle  corresponding  to  $ 


max 
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As  seen  from  Eq.  (113) ,  the  turn  rate  reaches  a  maximum 
at  an  angle  given  by 

=  cos-1 (P/2) 

with  a  corresponding  turn  rate  given  by 

V  (1+P/2)P  [1-(P/2)2]1_P/2 

ill  =  - 

1 v  1  max  K 

Thus,  the  value  of  K  =  K*  for  which  the  limiting  turn 
rate  is  achieved 

(1+P/2)P[1-(P/2)2]1_P/2 

K  =  S 

It  is  important  to  note  that,  if  ^o^ A'  the  "A" 

trajectory  is  invalid  since  the  maximum  turn  rate  of  the 
torpedo,  reached  at  the  range  r  =  3p  at  ♦  =<J>A  is  exceeded 
earlier  in  the  trajectory.  In  that  case,  the  limiting 
trajectory  is  the  trajectory.  Along  the  same  lines, 
we  argue  that  the  "B”  trajectory  always  is  a  limiting 
trajectory,  since  the  equality  (^0)B  cannot  be 

satisfied.  This  follows  from  the  fact  that 
0°«f>*  <  90° 

and  90°  as  4>B  is  the  supplement  of  <t>A- 

Below  some  critical  ship  speed,  the  torpedo  will  not 
be  turn  rate  limited.  This  speed  is  obtained  when  the 
trajectory  is  normal  to  the  ships  velocity  vector,  0=90°, 
at  a  range  of  r  ■  3p. 
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Again,  by  probing  Eq.  (88)  we  can  obtain  the  corres¬ 
ponding  values  of  <t>Q,  designated  (^Q)NL»  where  the  subscript 
"NL"  is  used  to  indicate  the  "no-limit"  boundary  point. 

The  computational  procedure  is  based  on: 

1.  Assuming  a  torpedo  speed  V  =  35  kts. 

2.  Using  an  initial  detection  range  (beginning  of 
homing)  rQ  =  250  m. 

3.  Assuming  the  following  maximum  turn  rates: 

a.  8  o/s  b.  12  o/s  c.  16  o/s 

d.  24  o/s  e.  36  o/s  f.  48  o/s 

4.  The  following  range  of  target  speeds 

0<V  <25  kts 
s 


The  aim  of  the  computation  is  to  determine  whether  the 
limiting  $  ,  and  hence  AOB,  is  governed  by  or  by  the 
turn  rate  at  the  range  r  -  3p. 


The  calculations  are  devided  into  two  parts,  and  are 
performed  on  a  TI  59  calculator. 

For  a  given  ship  speed,  Part  I  gives  the  sequential 
calculations  of  <t>A,  K^,  $B,  Kg,  0*  and  K*  for  each  of  the 
turn  rates.  The  computer  program  is  given  in  Appendix  E. 

Then,  in  Part  II,  the  probe  calculation  for  ( (J>Q ) A # 

(*  ) n,  (a  )*  =  ( a  n  *  are  performed.  The  program  is 
given  in  Appendix  E. 

The  limiting  results  are  given  on  a  polar  plot.  Fig.  41 

with 

<aob)a,b..  -  180°  -  ‘Va.b,* 

As  seen  from  Fig.  41,  for  a  target  speed  of  15  kts,  we 

need  AOB>85°  at  <b  for  a  maximum  turn  rate  16  o/s  in  order 

o 

to  have  a  hit  at  the  first  attack. 


IX.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  detection  performance  of  a  passive  homing  torpedo 
used  against  shallow-draft  surface  ships  operating  in 
moderate  sea  states  was  investigated.  Attention  was 
focused  on  the  effects  of  scattering  from  the  randomly 
rough  sea  surface  and  scattering  and  absorption  from  the 
bubble-dominated  inhomogeneous  layer  just  below  the  sea 
surface.  The  effects  of  these  two  scattering  mechanisms 
were  separately  estimated  and  their  relative  importance 
were  compared. 

The  passive  sonar  equation  was  used  to  predict  the 
performance  of  the  homing  system,  and  the  detection  range 
considering  these  two  scattering  effects  was  obtained  and 
compared  to  the  detection  range  based  on  a  reference 
model.  An  idealized  propagation  model  was  used  as 
reference  of  comparison.  This  reference  model  was  based 
on  a  noise  source  model  for  a  cavitating  propeller,  the 
operational  characteristics  for  a  square-law  detector, 
and  a  transmission  model  associated  with  a  homogeneous, 
unbounded  medium. 

Due  to  high  frequency,  moderate  sea  state  and  low 
grazing  angles,  the  scattering  from  the  randomly  rough 
sea  surface  was  found  to  be  dominated  by  the  direct  path. 
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This  result  includes  effects  from  geometrical 
shadowing . 

The  effect  of  scattering  and  absorption  from  the 
bubble-dominated,  inhomogeneous  subsurface  layer  was  in¬ 
vestigated  using  multiple  scattering  theory.  Both  the 
coherent  and  incoherent  limits  were  investigated  by  incor¬ 
porating  the  associated  absorption  coefficient  into  the 
transmission  equation. 

The  effect  on  the  sound  speed  from  the  bubble  content 
was  found  to  be  negligible. 

At  the  assumed  depth  setting  of  6  m  for  the  torpedo's 
search  and  attack  phase,  the  scattering  from  the  bubbles 
increased  the  transmission  loss.  This  increase  depended 
on  the  frequency  and  the  wind  speed. 

Two  operating  frequencies  were  investigated,  60  kHz 
and  30  kHz.  For  both  cases,  bubbles  significantly 
decreased  the  detection  range. 

For  a  torpedo  system  operating  at  the  high  frequencies, 
e.g. ,  60  kHz,  the  result  indicates  the  inadvisability  of 
using  a  searching  and  attack  depth  near  the  surface,  i.e., 
z  *  6  m. 

For  am  operating  frequency  around  30  kHz,  the  calculated 
detection  ramges  is  such  that  the  turn  rate  requirements 
for  a  pursuit  homing  trajectory  become  excessive.  For  a 
maximum  turn  rate  of  16  o/s,  this  limitation  can  be  avoided 


by  adapting  a  tactical  procedure  where  the  angle  on  bow 
(AOB)  at  the  beginning  of  the  torpedo  attack  is  greater 
than  85°. 

At  sea  state  3,  the  results  show  a  consistent  and 
general  trend  towards  the  need  for  lower  operating  fre¬ 
quency  in  order  to  increase  the  detection  range.  An 
operating  frequency  below  30  kHz  seems  indicated. 

Furthermore,  a  search  depth  below  the  bubble-dominated 
subsurface  layer,  i.e.,  z>15-20  m  would  result  in  an 
increased  detection  range. 

To  reduce  the  operational  limitations  induced  by  the 
scattering  and  absorption  effects,  a  high  maximum  turn 
rate  together  with  a  variable  speed  capability,  where 
p^2  would  be  beneficial. 

The  result  of  this  analysis  has  clearly  demonstrated 
the  importance  of  environmental  factors  on  the  torpedo 
capability,  and  is  useful  in  giving  insight  into  the  be¬ 
havior  of  a  homing  torpedo  during  its  critical  attack 
phase . 

A  valuable  follow-on  of  this  study  would  be  an  investi¬ 
gation  of  the  effects  of  the  bubble-dominated  subsurface 
layer  on  target  validation  and  pitch  plane  steering  when 
the  torpedo  search  depth  is  50-60  m. 
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APPENDIX  A 


9 

DETAILED  OCEANOGRAPHIC  BACKGROUND  MATERIAL 

A.  GEOGRAPHIC  DESCRIPTION 
The  Norwegian  coastal  waters  constitute  the  eastern 

boundary  of  the  Norwegian  Sea.  Although  some  general 
aspects  related  to  the  Norwegian  Sea  will  be  covered, 
this  analysis  will  be  concentrated  on  the  Norwegian 
coastal  waters  above  68°N. 

B.  WIND 

The  northern  region  of  the  Norwegian  Sea  is  affected 
by  the  Polar  Easterlies  and  the  southern  region  by  the 
prevailing  Westerlies.  There  are  two  dominant  air  masses 
which  are  relatively  permanent: 

-The  Greenland  high. 

-The  Iceland  low. 

These  pressure  systems  produce  storms  which  are  carried 
across  the  Norwegian  Sea  in  a  belt  from  Iceland  towards 
the  Norwegian  Coast  causing  steady  precipitation  and  wind 
most  of  the  year.  The  steep  Norwegian  Coast  has  a  con¬ 
siderable  influence  on  the  winds  and  consequently  also  on 
the  waves  in  the  coastal  waters.  The  main  general  modi¬ 
fications  are  that  the  streamlines  tend  to  run  parallel 
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to  the  coast  and  that  wind  and  sea  increases  with  distance 
from  the  sheltered  coast  into  open  ocean. 

Strong  local  variation  may  occur.  The  most  important 
of  these  are  the  marked  local  increase  in  wind  speed  in 
areas  where  the  coast  sharply  changes  direction.  One  such 
"corner  effect,"  caused  by  the  confluence  of  the  streamlines, 
occurs  near  "Nordkapp"  (North  Cape) .  Also  of  importance 
are  monsoonal  effects  due  to  the  different  heat  capacities 
between  the  ocean  and  the  continent.  Drainage  of  cold  air 
from  the  inland  valleys  in  the  wintertime  causes  a  marked 
increase  in  the  wind  speed  in  several  areas  along  the 
coast.  Most  of  these  coastal  effects  are  significantly 
dissipated  at  distances  of  approximately  50  nmi.  from  the 
coast.  The  fact  that  the  wind  tends  to  blow  along  the 
coast  is  clearly  demonstrated  in  Figs.  12  and  13  which 
include  data  from  weather  stations  from  "Hillesoy"  to 
"Ona"  and  "Myken"  to  "Furuholmen"  respectively.  The  high 
frequency  of  offshore  winds  is  caused  by  the  drainage  of 
cold  air  from  the  inland  valleys  during  winter  time. 

That  this  phenomenon  is  closely  connected  to  the  coast 
is  illustrated  by  the  fact  that  it  is  missing  at  weather 
station  "Skomvaer"  situated  approximately  50  nmi.  off  the 
main  coast.  A  frequency  distribution  of  observed  wind 
speeds  along  the  coast,  obtained  from  Ref.  1,  is  presented 
in  Table  V.  A  summary  of  this  table  follows: 
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TABLE  V 


FREQUENCY  DISTRIBUTION  OF 

WIND 

SPEED  IN 

PERCENT 

PER  YEAR 

AT  WEATHER  STATIONS  ALONG 

THE 

NORWEGIAN 

COAST 

m/sec 

0-7 

8-13 

14-20 

*21 

Beaufort 

0-4 

5+6 

7+8 

=  9 

Ferder 

67 

28 

5 

0.1 

Lyngor 

79 

19 

2 

0.2 

Lista 

67 

28 

5 

0.2 

Utsira 

73 

22 

5 

0.3 

Hellisoy 

77 

19 

4 

0.2 

Krakenes 

58 

28 

11 

2.7  Jan-Dec 

1949-1975 

Ona 

70 

22 

7 

u  *  / 

Sula 

61 

30 

8 

1.1 

Nordoyan 

49 

35 

14 

1.9 

Myken 

64 

26 

9 

1.2 

Skomvaer 

58 

31 

10 

0.9 

Andenes 

79 

18 

3 

0.1 

Torsvar 

72 

23 

5 

0.4 

Fruholmen 

54 

32 

12 

1.8 

Vardo 

75 

22 

3 

0.1 

Biomoya 

63 

31 

6 

0.4  1956-1975 

Polarfront 

46 

30 

14 

1.4 
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The  highest  winds  are  reported  from  the  areas  between 
62°N  and  68°N.  The  frequency  of  high  winds  in  this  area 
are  significantly  greater  than  those  found  at  stations  to 
the  south  of  62°N  and  also  to  the  north  of  68°N.  The  wind 
condition  around  North  Cape,  are  very  severe.  This  area 
can  be  compared  with  the  other  coastal  area  of  high  wind 
speed  such  as  "Stadt"  (represented  by  the  weather  station, 
"Krakenes") .  In  these  areas  the  frequency  of  storms  is 
greater  than  at  "Polar front"  situated  in  the  open  Norwegian 
Sea. 


C .  WAVES 

Frequency  distribution  of  significant  wave  heights  are 
represented  in  Table  VI. 

The  station  north  of  68°N  is  characterized  by  compara¬ 
tively  small  frequency  of  high  waves.  Even  at  "Furuholmen” 
where  wind  conditions  are  very  severe,  the  frequency  dis¬ 
tribution  of  significant  wave  heights  is  similar  to  more 
sheltered  areas  like  "Utsira." 

The  seasonal  variation  for  the  area  of  interest  around 
70°N  is  given  in  Figs.  3  and  4.  The  average  monthly  dis¬ 
tribution  of  significant  wave  heights  for  a  typical 
station  like  "Andenes"  is  given  in  Fig.  3  and  Fig.  4 
gives  the  yearly  distribution  for  this  station. 
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TABLE  VI 


FREQUENCY  DISTRIBUTION  OF  SIGNIFICANT  WAVE 
HEIGHT  IN  PERCENT  AT  THE  WEATHER  STATION  "ANDENES" 


STATION:  Andenes  YEARS:  1949-1972 


SS 

0+1+2 

3 

4 

5 

6 

7 

8 

9 

Hc 

0 

0.5  - 

1.3 

-2.5- 

4.0 

6.0 

-  9.0- 

13.9 

N 

s 

0.4 

1.2 

2.4 

3.9 

5.9 

8.9 

13.9 

JAN 

15.3 

26.9 

33.3 

15.5 

6.6 

2.3 

0.18 

0.00 

2542 

FEB 

17.7 

28.7 

31.3 

14.2 

5.9 

2.0 

0.14 

0.00 

2295 

MAR 

18.1 

28.9 

30.8 

14.3 

5.8 

1.9 

0.17 

0.00 

2529 

APR 

24.5 

33.1 

27.4 

10.8 

3.4 

0.8 

0.04 

2293 

MAY 

30.9 

35.5 

23.9 

7.7 

1.7 

0.3 

0.01 

1720 

JUN 

37.4 

36.9 

20.2 

4.8 

0.6 

0.1 

0.00 

1565 

JUL 

42.6 

36.2 

17.1 

3.7 

0.4 

0.0 

0.00 

1618 

AUG 

40.6 

36.4 

18.0 

4.3 

0.5 

0.1 

1621 

SEP 

29.7 

34.3 

24.4 

8.9 

2.1 

0.5 

0.01 

2301 

OCT 

21.8 

31.3 

29.4 

12.4 

4.1 

1.0 

0.01 

2374 

NOV 

21.7 

31.5 

29.9 

11.7 

4.1 

1.1 

0.06 

2303 

DEC 

18.5 

30.1 

32.2 

13.2 

4.7 

1.2 

0.09 

2379 

YEAR 

26.6 

32.5 

26.5 

10.1 

3.3 

0.9 

0.06 

0.00 

N  =  Number  of  Observations 
SS  =  Class  Interval,  State  of  Sea 
Hg  =  Significant  Wave  Height 
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D.  AMBIENT  NOISE 


Due  to  the  fact  that  the  Norwegian  Sea  is  physically 
separated  from  the  Atlantic  by  the  Faeroy-Shetland-Iceland 
ridge,  and  from  the  Greenland  Sea  by  the  Jan  Mayan  ridge, 
little  long  distance  shipping  noise  is  transferred  into 
the  area.  This,  combined  with  relatively  low  shipping 
traffic  in  the  central  and  northern  parts  of  the  Norwegian 
Sea,  produces  a  relatively  low  ambient  noise  level  for  the 
frequency  band  100<f<1000  Hz,  especially  when  noise  from 
marine  life  is  not  included.  Further,  in  the  absence  of 
nearby  shipping  and  marine  life,  the  ambient  noise  level 
in  the  frequency  band  l<f<50  kHz  is,  according  to  Ref.  4, 
dominated  by  the  wind. 

A  typical  area  for  the  central  part  of  the  Norwegian 
Sea  can  be  represented  by  the  weather  station  "Polarfront" 
at  66°N,  2°E.  Reference  1  shows  that  there  is  approximately 
a  15%  chance  of  finding  wind  forces  of  Beaufort>6.  These 
effects  predict  a  moderate  ambient  noise  level  in  the 
frequency  range  100  Hz-50  kHz  for  the  central  and  northern 
regions  of  the  Norwegian  Sea. 

The  shallow  coastal  Norwegian  waters  are,  according  to 
Ref.  4,  typically  5-10  dB  noisier  than  the  corresponding 
deep  waters.  However,  great  variability  caused  by  local 
traffic,  fishing  fleet  activity,  marine  life,  and  local 
wind  conditions  makes  ambient  noise  level  prediction 
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difficult  in  these  areas.  This  means  that  accurate  ambient 
noise  level  determinations  have  to  be  made  on  the  spot,  as 
it  is  both  site  and  time  dependent. 

A  noticeable  influence  on  the  ambient  noise  level  is 
rain,  which  is  a  year  around  feature  along  the  Norwegian 
coast.  As  seen  from  Fig.  8,  taken  from  Ref.  4,  rain  has 
a  tendency  to  produce  a  constant  high  ambient  noise  level 
over  a  large  frequency  range,  thus  dominating  ot7  -?r 
effects.  Furthermore,  for  the  upper  frequency  of  interest, 
i.e.,  around  60  kHz  the  lower  bound  for  the  ambient 
noise  is  determined  by  the  thermal  agitation,  see  Fig.  8. 

In  determining  the  figure  of  merit  (FOM)  for  a 
passive  sonar  system,  the  noise  level  will  be  the  larger 
of  either  the  self  noise  or  the  ambient  noise.  For  a 
torpedo  the  self  noise  will  typically  be  dominant. 

E.  SOUND  SPEED  PROFILES 

Again,  concentrate  on  data  relevant  to  Norwegian 
coastal  waters.  According  to  Ref.  6,  which  covers  the 
southern  part  of  the  Norwegian  coast,  low  sound  speeds 
are  common  because  of  the  influence  of  water  from  the 
Baltic  Sea  combined  with  fresh  water  drainage  from  the 
fjords.  Furthermore,  great  variability,  both  seasonal 
and  within  seasons,  is  encountered.  Figure  9  obtained 
from  Ref.  5  gives  a  picture  of  the  sound  speed  profiles 
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for  the  northern  Norwegian  coast.  Again,  large  variations 
are  common.  Noticeable  in  both  sets  of  data  is  a  typical 
seasonal  pattern  of  strong  cooling  of  the  surface  layer 
during  winter  and  a  similarly  strong  heating  during 
summer.  Furthermore,  note  that  the  minimum  and  maximum 
are  relatively  shallow,  i.e.,  less  than  50  m. 

Also  characteristic  is  the  influence  of  the  cold  and 
fresh  melt  waters  drained  out  through  the  fjord-arms 
during  spring  and  summer. 

To  illustrate  the  sonar  problems  associated  with  these 
sound  speed  profiles,  ray  paths  for  the  extremes  of  Fig.  9 
are  shown  in  Figs.  10  and  11,  where  the  source  is  3  m  below 
the  sea  surface. 
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APPENDIX  B 


SURFACE  SCATTERING  TI  59  PROGRAM 

A.  INTRODUCTION 

This  program  gives  specular  scattered  power  at  the 
receiver  versus  incoming  power  at  the  randomly  rough 
surface  in  the  high  frequency  limit  according  to  Eq.  (70) . 

Shadowing  of  surface  areas  by  other  parts  of  the 
boundary  are  taken  into  account  by  the  bistatic  shadowing 
function  S(0). 

Furthermore,  the  program  gives  the  effect  of  the 
randomly  rough  surface  compared  to  the  idealized  free- 
field  condition  as  expressed  in  Eg.  (80).  The  results 
of  these  calculations  is  given  in  Fig.  27. 

B .  PROGRAM  STEPS 

A  block-diagram  of  the  computer  program  is  given  in 
Fig.  26. 

The  program  uses  the  partitioning  ratio  of  program 
to  data  space  according  to  code  4  OP17.  The  users  in¬ 
structions  are  as  follows: 
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Procedure 

Enter 

Press 

Display 

Enter  data 

Detection  range 

2nd  A 

R  1 

rd 

Enter  data 

Source  depth 

R/S 

h  1 
s 

Enter  data 

Receiver  depth 

R/S 

h  1 
r 

Enter  data 

Beam  width 

R/S 

(f)1 

Enter  data 

Wind  speed 

R/S 

w1 

Calculate  0 , 

Ri ,R2 

— 

2nd  B 

A 

( 0  ,  Ri ,  R2  ,  A) 

Calculate  <i;12> 

— 

2nd  C 

<^12>1 

Calculate  v 

— 

2nd  D 

v1 

Calculate  erfc  v 

— 

A 

erfc  v1 

Calculate  S(0) 

— 

B 

SO)1 

<s2  > 

Calculate  <ptPj*> 

— 

C 

<s2>1 

^iPi^ 

Calculate  AIL 

— 

D 

AIL1 

■'"These  values  are  printed  automatically  if  the  calculator 
is  connected  to  the  PC-100A  Print  Cradle. 


For  the  error  function  complement  we  have 


erfc(v)  =  —  f  e  a  da 
/?  v 

in  contrast  to  the  normal  distribution 


Q(u)  =*  —  /  e"fc2/2  dt 

m  a 

However,  there  is  a  linear  relationship  between  the 


two  functions 
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The  numerical  equation  used  to  calculate  the  erfc(v) 
is  a  modified  program  from  Texas  Instruments  [Ref.  24]. 
erfc(v)  =  Z  (v)  [bia+b2Ct2+b  jet  3+bi*ot‘*+bsa5  ] 

where 

a  =  I+pv 

p  =  .231649 

b  =  .451673691 
i 

b2  =  -.504257336 
b3  =  2.51939026 
b4  =  -2.563346623 
bs  =  1.881292139 


The  program  steps  are  listed  below,  giving  location 


(LOC) ,  code  (COD) ,  key  symbol  (KEY) ,  and  comments. 


000 

76 

LBL 

023 

17 

B' 

046 

05 

05 

001 

16 

A' 

024 

53 

( 

047 

55 

T 

002 

42 

ST 

025 

53 

( 

048 

43 

RCL 

003 

00 

00 

026 

43 

RCL 

049 

02 

02 

004 

99 

PRT 

027 

00 

00 

050 

54 

) 

005 

91 

R/S 

028 

65 

X 

051 

22 

INV 

006 

42 

STO 

029 

43 

RCL 

052 

30 

TAN 

007 

01 

01 

030 

02 

02 

053 

54 

) 

008 

99 

PRT 

031 

54 

) 

054 

42 

STO 

009 

91 

R/S 

032 

55 

+ 

055 

06 

06 

010 

42 

STO 

033 

53 

( 

056 

99 

PRT 

Oil 

02 

02 

034 

43 

RCL 

057 

53 

( 

012 

99 

PRT 

035 

01 

01 

058 

53 

( 

013 

91 

R/S 

036 

85 

+ 

059 

43 

RCL 

014 

42 

STO 

037 

43 

RCL 

060 

05 

05 

015 

03 

03 

038 

02 

02 

061 

55 

T 

016 

99 

PRT 

039 

54 

) 

062 

53 

( 

017 

91 

R/S 

040 

54 

) 

063 

43 

RCL 

018 

42 

STO 

041 

42 

STO 

064 

06 

06 

019 

04 

04 

042 

05 

05 

065 

38 

SIN 

020 

99 

PRT 

043 

53 

( 

066 

54 

) 

021 

91 

R/S 

044 

53 

( 

067 

54 

) 

022 

76 

LBL 

045 

43 

RCL 

068 

42 

STO 

112 


I 

r 


r 


i 

i 


I 


069 

07 

07 

120 

43 

RCL 

171 

69 

OP 

070 

99 

PRT 

121 

03 

03 

172 

06 

06 

071 

53 

( 

122 

38 

SIN 

173 

91 

R/S 

072 

53 

( 

123 

65 

X 

174 

76 

LBL 

073 

43 

RCL 

124 

43 

RCL 

175 

19 

D' 

074 

00 

00 

125 

07 

07 

176 

53 

( 

075 

75 

- 

126 

54 

) 

177 

43 

RCL 

076 

43 

RCL 

127 

54 

) 

178 

06 

06 

077 

05 

05 

128 

42 

STO 

179 

30 

TAN 

078 

54 

) 

129 

10 

10 

180 

35 

1/x 

079 

55 

T 

130 

69 

OP 

181 

54 

) 

080 

53 

( 

131 

06 

06 

182 

42 

STO 

081 

43 

RCL 

132 

91 

R/S 

183 

12 

12 

082 

06 

06 

133 

76 

LBL 

184 

53 

( 

083 

38 

SIN 

134 

18 

C' 

185 

43 

RCL 

084 

54 

) 

135 

25 

CLR 

186 

12 

12 

085 

54 

) 

136 

69 

OP 

187 

55 

T 

086 

42 

STO 

137 

00 

00 

188 

53 

( 

087 

08 

08 

138 

00 

00 

189 

53 

( 

088 

99 

PRT 

139 

06 

6 

190 

02 

2 

089 

53 

( 

140 

02 

2 

191 

65 

X 

090 

43 

RCL 

141 

07 

7 

192 

43 

RCL 

091 

07 

07 

142 

03 

3 

193 

11 

11 

092 

85 

+ 

143 

02 

2 

194 

54 

) 

093 

43 

RCL 

144 

03 

3 

195 

34 

>6c 

094 

08 

08 

145 

03 

3 

196 

54 

) 

095 

54 

) 

146 

69 

OP 

197 

54 

) 

096 

42 

STO 

147 

04 

04 

198 

42 

STO 

097 

09 

09 

148 

53 

( 

199 

13 

13 

098 

25 

CLR 

149 

53 

( 

200 

99 

PRT 

099 

69 

OP 

150 

53 

( 

201 

91 

R/S 

100 

00 

00 

151 

43 

RCL  . 

202 

76 

LBL 

101 

01 

1 

152 

04 

04 

203 

11 

A 

102 

03 

3 

153 

65 

X 

204 

53 

( 

103 

03 

3 

154 

05 

5 

205 

53 

( 

104 

05 

5 

155 

93 

■ 

206 

53 

( 

105 

01 

1 

156 

01 

1 

207 

53 

( 

106 

07 

7 

157 

02 

2 

208 

43 

RCL 

107 

01 

1 

158 

54 

) 

209 

13 

13 

108 

03 

3 

159 

85 

+ 

210 

33 

X2 

109 

69 

OP 

160 

03 

3 

211 

54 

( 

110 

04 

04 

161 

54 

) 

212 

55 

T 

111 

53 

( 

162 

65 

X 

213 

02 

2 

112 

53 

( 

163 

93 

. 

214 

54 

) 

113 

89 

IT 

164 

00 

0 

215 

94 

+/- 

114 

65 

X 

165 

00 

0 

216 

22 

INV 

115 

43 

RCL 

166 

01 

1 

217 

23 

LNX 

116 

08 

08 

167 

54 

) 

218 

54 

) 

117 

54 

) 

168 

50 

lxl 

219 

65 

X 

118 

65 

X 

169 

42 

STO 

220 

53 

( 

119 

53 

( 

170 

11 

11 

221 

53 

( 

;  t 


i  : 

-  i 


i 

* 

i 


113 


•f  »1»"*  «f 


T 


I 

! 


222  89  -  273 

223  34  /x  274 

224  65  x  275 

225  02  2  276 

226  54  )  277 

227  35  1/X  278 

228  54  )  279 

229  54  )  280 

230  42  STO  281 

231  14  14  282 

232  53  (  283 

!  233  53  (  284 

j  234  53  (  285 

r  235  53  (  286 

236  43  RCL  287 

237  13  13  288 

j  238  65  X  289 

I  239  93  .  290 

240  02  2  291 

241  03  3  292 

242  01  1  293 

243  06  6  294 

244  04  4  295 

!  245  01  1  296 

246  09  9  297 

247  86  +  298 

248  01  1  299 

249  54  )  300 

250  35  1/ x  301 

251  54  )  302 

252  42  STO  303 

253  15  15  304 

254  53  (  305 

i  255  53  (  306 

256  43  RCL  307 

257  14  14  308 

258  65  x  309 

259  01  1  310 

260  93  .  311 

261  08  8  312 

262  01  1  313 

263  54  )  314 

264  65  x  315 

265  53  {  316 

266  43  RCL  317 

267  15  15  318 

268  65  x  319 

269  93  .  320 

270  03  3  321 

271  .01  1  322 

272  09  9  323 


I 

ii 


03  3 

324 

43  RCL 

08  8 

325 

15  15 

01  1 

326 

45  Yx 

05  5 

327 

04  4 

03  3 

328 

85  + 

75  - 

329 

01  1 

93  . 

330 

93  • 

03  3 

331 

03  3 

05  5 

332 

03  3 

06  6 

333 

00  0 

05  5 

334 

02  2 

06  6 

335 

07  7 

03  3 

336 

04  4 

07  7 

337 

04  4 

08  8 

338 

02  2 

02  2 

339 

09  9 

65  x 

340 

00  0 

43  RCL 

341 

65  X 

15  15 

342 

43  RCL 

33  x2 

343 

15  15 

85  + 

344 

45  Yx 

01  1 

345 

05  5 

93  . 

346 

54  ) 

97  7 

347 

54  ) 

08  8 

348 

42  STO 

01  1 

349 

16  16 

04  4 

350 

99  PRT 

07  7 

351 

91  R/S 

07  7 

352 

76  LBL 

09  9 

353 

12  B 

03  3 

354 

53  ( 

07  7 

355 

53  ( 

65  x 

356 

53  ( 

43  RCL 

357 

43  RCL 

15  15 

358 

13  13 

45  Yx 

359 

33  x2 

03  3 

360 

94  +/- 

75  - 

361 

22  INV 

01  1 

362 

23  LNX 

93  . 

363 

54  ) 

08  8 

364 

75  - 

02  2 

365 

53  ( 

01  1 

366 

89  n 

02  2 

367 

34  /x 

05  5 

368 

65  x 

05  5 

369 

43  RCL 

09  9 

370 

13  13 

07  7 

371 

65  x 

08  8 

372 

43  RCL 

00  0 

373 

16  16 

65  x 

374 

54  ) 

114 

APPENDIX  C 


BUBBLE  DYNAMICS  TI  59  PROGRAM 

A.  INTRODUCTION 

This  program  gives  the  resonant  frequency  f  and  the 
damping  constant  <5  for  bubbles  according  to  Eq.  (92) . 
Furthermore,  the  program  gives  scattering  cross  section 
a  s,  extinction  cross  section  oe,  and  absorption  cross 
section  a  as  a  function  of  bubble  radius  a,  incoming  fre- 

a 

quency  f,  and  depth  z.  The  results  of  the  calculations 
are  given  in  Figs.  31  and  32. 


B.  PROGRAM  STEPS 


The  user ' s 

instructions  are  as 

follows: 

Procedure 

Enter 

Press 

Display 

Enter  data 

Bubble  Radius 

2nd  A 

a1 

Enter  data 

Incoming  Frequency 

R/S 

f1 

Enter  data 

Depth 

R/S 

z1 

Calculate  X 
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2nd  B 
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Calculate 
coshX  and 
sinhX 

2nd  C 
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sinhX1 
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— 

2nd  D 

d/b1 

Calculate  b 
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2nd  E 

b1 

Calculate  8 

— 
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B1 

Calculate  f 

— 
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V 

Calculate  6 

— 

C 

51 

Calculate  a 

— 

D 

ae1 

Calculate  a 

and  o  8 

s 

— 

E 

Os1 

“a1 
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^These  values  are  printed  automatically  if  the  calculator 
is  connected  to  PC-100A  Print  Cradle. 

The  program  steps  are  listed  below  giving  location  (LOC) , 


code  (COD) ,  key  (KEY)  and  comments. 
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43 
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53 
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APPENDIX  D 


NUMERICAL  INTEGRATION  TI  59  PROGRAM 


A-  INTRODUCTION 

In  order  to  perform  the  numerical  integration  of  an (a) da, 
a  standard  Texas  Instrument's  program  was  used  [Ref.  25]. 

This  program  performs  the  integration  by  using  Simpson's 
discrete  approximation  based  on  the  following  expression 
x  . 

I  =  f  n  f (x) dx  =  3 (fQ+4f i+2f 2+4f 3+2f 4  + _ +2 f n_  +4 f  a  +f n > 


where 


f(x)  must  be  known  at  n+1  equally  spaced  points  (fQ-fn) . 

x  -x_ 
h  =  -2 — ° 
n 


'  n  o 


n+7  -  number  of  data  registers  available 
n = number  of  subintervals  =2,  4,  6 


B.  PROGRAM  INSTRUCTIONS 

The  program  is  taken  from  the  master  library  program 
package  by  using  the  code  2nd  Pgm  10  on  the  calculator. 
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The  user’s  instructions  are  as  follows 


Procedure 

Enter 

Press 

Display 

Enter  data 

Subintervals 

A 

n1 

Enter  data 

h 

B 

h1 

Enter  data 

Function  values 

0 

C 

0 

fo 

R/S 

fl 

R/S 

• 

V 

fn 

• 

• 

R/S 

f  1 
n 

Calculate 

- 

D 

I1 

■''These  values  are  printed  automatically  if  the  calculator 
is  connected  to  PC-100A  Print  Cradle. 


The  program  steps  are  listed  below  giving  location  (LOC) , 


key  symbol  (KEY) ,  and  comments. 
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01 
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004 
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023 
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32 
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005 
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35 
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98 
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55 
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044 

92 

RTN 

007 

02 
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LBL 

045 

76 

LBL 

008 

54 

) 

027 

12 

B 

046 

50 

lxl 

009 

42 

STO 

028 

42 

STO 

047 

76 

ST* 

010 

02 

02 

029 

03 

03 

048 

01 

01 

011 

22 

INV 

030 

99 

PRT 

049 

32 

x£t 

012 

59 

INT 

031 

92 

RTN 

050 

01 

1 

013 

29 

CP 

032 

76 

LBL 

051 

44 

SUM 

014 

22 

INV 

033 

13 

C 

052 

01 

01 

015 

67 

EQ 

034 

53 

( 

053 

32 

x$t 

016 

52 

EE 

035 

24 

CE 

054 

99 

PRT 

017 

43 

RCL 

036 

85 

+ 

055 

92 

RTN 

018 

05 

05 

037 

32 

x  >  t 

056 

61 

GTO 
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APPENDIX  E 


TURN  RATE  LIMITATION  TI  59  PROGRAMS 


A.  INTRODUCTION 

The  turn  rate  limitation  calculations  are  divided  into 
two  parts  with  separate  programs. 

B.  PART  I  PROGRAM 

For  a  given  initial  detection  range  and  a  given  ship 
speed.  Part  I  performs  a  sequential  calculation  of 
<(>  ,  ko  ,  (J>*  and  K*  for  each  of  the  maximum  turn  rates  investi- 
gated,  together  with  the  "no  limit"  conditions  (VS)NL  and 
based  on 


K* 


V  (l+P/2)P[l-(P/2) 2] 1"P/2 
s 


<Vnl  * 

knl  ■  3p 
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For  p>2,  the  solution  for  *  and  K*  are  not  valid. 
Furthermore,  for  Vg  (VS)NL  no  solutions  are  valid  for  any 
of  the  quantities. 

The  program  is  based  on  the  fixed  torpedo  speed  of 
VT  =  35  kts  (18  m/s)  and  an  initial  detection  range  of 

To  3  -  "d  *  250  *>• 

The  user  instructions  for  the  program  are  as  follows: 


Prodcedure 

Enter 

Press 

Display 

Enter  data 

Maximum  Turn  Rate 

A 

U|  1 

1 Y 'max 

Enter  data 

Ship  speed 

B 

V  1 

T 

Calculate  <$* 
and  K* 

— 

C 

Calculate  <J> 
and  4>b 

— 

D 

<fr  ^R1 

A'  B 

Calculate  K 
and  Kb  a 

—  — 

E 

Vb1 

^These  values  are  printed  automatically  if  the  calculator 
is  connected  to  PC-100A  Print  Cradle. 


A  listing  of  the  program  steps  follows  giving  location 


(LOC) ,  code  (COD),  key  symbol  (KEY),  and  comments. 


000 

76 

LBL 

012 

69 

OP 

024 

54 

) 

001 

11 

A 

013 

00 

00 

025 

42 

STO 

002 

42 

STO 

014 

03 

3 

026 

03 

03 

003 

01 

01 

015 

03 

3 

027 

69 

OP 

004 

99 

PRT 

016 

69 

OP 

028 

06 

06 

005 

91 

R/S 

017 

04 

04 

029 

25 

CLR 

006 

76 

LBL 

018 

53 

( 

030 

69 

OP 

007 

12 

B 

019 

01 

1 

031 

00 

00 

008 

42 

STO 

020 

08 

8 

032 

02 

2 

009 

02 

02 

021 

55 

f 

033 

06 

6 

010 

99 

PRT 

022 

43 

RCL 

034 

03 

3 

Oil 

25 

CLR 

023 

02 

02 

035 

01 

1 

125 


036 

02 

2 

087 

03 

3 

138 

45  ■ 

yx 

037 

07 

7 

088 

03 

3 

139 

43 

RCL 

038 

69 

OP 

089 

02 

2 

140 

03 

03 

039 

04 

04 

090 

03 

3 

141 

54 

) 

040 

53 

( 

091 

05 

5 

142 

54 

) 

041 

43 

RCL 

092 

01 

1 

143 

65 

X 

042 

03 

03 

093 

69 

OP 

144 

53 

( 

043 

65 

X 

094 

04 

04 

145 

53 

( 

044 

03 

3 

095 

53 

( 

146 

01 

1 

045 

54 

) 

096 

53 

( 

147 

75 

— 

046 

42 

STO 

097 

53 

( 

148 

53 

( 

047 

10 

10 

098 

43 

RCL 

149 

53 

( 

048 

69 

OP 

099 

03 

03 

150 

43 

RCL 

049 

06 

06 

100 

55 

T 

151 

03 

03 

050 

25 

CLR 

101 

02 

2 

152 

55 

T 

051 

69 

OP 

102 

54 

) 

153 

02 

2 

052 

00 

00 

103 

22 

INV 

154 

54 

) 

053 

04 

4 

104 

39 

COS 

155 

33 

X2 

054 

02 

2 

105 

54 

) 

156 

54 

) 

055 

03 

3 

106 

54 

) 

157 

54 

) 

056 

06 

6 

107 

42 

STO 

158 

45 

yX 

057 

03 

3 

108 

04 

04 

159 

53 

( 

058 

01 

1 

109 

69 

OP 

160 

01 

1 

059 

02 

2 

110 

06 

06 

161 

75 

— 

060 

07 

7 

111 

25 

CLR 

162 

53 

( 

061 

69 

OP 

112 

69 

OP 

163 

43 

RCL 

062 

04 

04 

113 

00 

00 

164 

03 

03 

063 

53 

( 

114 

02 

2 

165 

55 

T 

064 

53 

( 

115 

06 

6 

166 

02 

2 

065 

53 

( 

116 

05 

5 

167 

54 

) 

066 

03 

3 

117 

01 

1 

168 

54 

) 

067 

65 

X 

118 

69 

OP 

169 

54 

) 

068 

01 

1 

119 

04 

04 

170 

54 

) 

069 

08 

8 

120 

53 

( 

171 

55 

T 

070 

54 

) 

\21 

53 

( 

172 

43 

RCL 

071 

65 

X 

122 

53 

( 

173 

01 

01 

072 

43 

RCL 

123 

53 

( 

174 

54 

) 

073 

01 

01 

124 

43 

RCL 

175 

42 

STO 

074 

54 

) 

125 

02 

02 

176 

05 

05 

075 

34 

/£ 

126 

65 

X 

177 

69 

OF 

076 

54 

) 

127 

53 

( 

178 

06 

06 

077 

42 

STO 

128 

53 

( 

179 

91 

R/S 

078 

11 

11 

129 

01 

1 

180 

76 

LBL 

079 

69 

OP 

130 

85 

+ 

181 

14 

D 

080 

06 

06 

131 

53 

( 

182 

25 

CLR 

081 

91 

R/S 

132 

43 

RCL 

183 

69 

OP 

082 

76 

LBL 

133 

03 

03 

184 

00 

00 

083 

13 

C 

134 

55 

T 

185 

03 

3 

084 

25 

CLR 

135 

02 

2 

186 

03 

3 

085 

69 

OP 

136 

54 

) 

187 

02 

2 

086 

00 

00 

137 

54 

) 

188 

03 

3 

|4f> d  *'  »"»' 


1 


189 

01 

1 

240 

06 

06 

291 

54 

) 

190 

03 

3 

241 

91 

R/S 

292 

54 

) 

191 

69  < 

OP 

242 

76 

LBL 

293 

54 

) 

192 

04 

04 

243 

15 

E 

294 

42 

STO 

193 

53 

( 

244 

25 

CLR 

295 

08 

08 

194 

53 

( 

245 

69 

OP 

296 

69 

OP 

195 

53 

( 

246 

00 

00 

297 

06 

06 

196 

53 

( 

247 

02 

2 

298 

25 

CLR 

197 

53 

( 

248 

06 

6 

299 

69 

OP 

198 

43 

RCL 

249 

01 

1 

300 

00 

00 

199 

03 

03 

250 

03 

3 

301 

02 

2 

200 

65 

X 

251 

69 

OP 

302 

06 

6 

201 

03 

3 

252 

04 

04 

303 

01 

1 

202 

54 

) 

253 

53 

( 

304 

04 

4 

203 

65 

X 

254 

53 

( 

305 

69 

OP 

204 

43 

RCL 

255 

43 

RCL 

306 

04 

04 

205 

01 

01 

256 

02 

02 

307 

53 

( 

206 

54 

) 

257 

65 

X 

308 

53 

( 

207 

55 

258 

53 

( 

309 

43 

RCL 

208 

43 

RCL 

259 

53 

( 

310 

02 

02 

209 

02 

02 

260 

01 

1 

311 

65 

X 

210 

54 

) 

261 

85 

+ 

312 

53 

( 

211 

22 

INV 

262 

53 

( 

313 

53 

( 

212 

38 

SIN 

263 

43 

RCL 

314 

01 

1 

213 

54 

) 

264 

06 

06 

315 

85 

+ 

214 

42 

STO 

265 

39 

COS 

316 

53 

( 

215 

06 

06 

266 

54 

> 

317 

43 

RCL 

216 

69 

OP 

267 

54 

) 

318 

07 

07 

217 

06 

06 

268 

45 

yX 

319 

39 

COS 

218 

25 

CLR 

269 

43 

RCL 

320 

54 

) 

219 

69 

OP 

270 

03 

03 

321 

54 

>x 

220 

00 

00 

271 

54 

) 

322 

45 

y 

221 

03 

3 

272 

54 

) 

323 

43 

RCL 

222 

03 

3 

273 

55 

T 

324 

03 

03 

223 

02 

2 

274 

53 

( 

325 

54 

) 

224 

03 

3 

275 

43 

RCL 

326 

54 

) 

225 

01 

1 

276 

01 

01 

327 

55 

T 

226 

04 

4 

277 

65 

X 

328 

53 

( 

227 

69 

OP 

278 

53 

( 

329 

43 

RCL 

228 

04 

04 

279 

53 

( 

330 

01 

01 

229 

53 

( 

280 

43 

RCL 

331 

65 

X 

230 

01 

1 

281 

06 

06 

332 

53 

( 

231 

08 

8 

282 

38 

SIN 

333 

53 

( 

232 

00 

0 

283 

54 

>x 

334 

43 

RCL 

233 

75 

284 

45 

Y 

335 

07 

07 

234 

43 

RCL 

285 

53 

( 

336 

38 

SIN 

235 

06 

06 

286 

43 

RCL 

337 

54 

^  X 

236 

54 

) 

287 

03 

03 

338 

45 

y 

237 

42 

STO 

288 

75 

- 

339 

53 

( 

238 

07 

07 

289 

02 

2 

340 

43 

RCL 

239 

69 

i  OP 

290 

54 

i  ) 

341 

03 

03 

rr 

j. 
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342 

75 

- 

343 

02 

2 

344 

54 

) 

345 

54 

) 

346 

54 

) 

347 

54 

) 

348 

42 

STO 

349 

09 

09 

350 

69 

OP 

351 

06 

06 

352 

91 

R/S 

353 

02 

2 

354 

54 

) 

355 

54 

) 

356 

54 

) 

357 

54 

) 

358 

42 

STO 

359 

09 

09 

360 

69 

OP 

361 

06 

06 

362 

91 

R/S 

C.  PART  II  PROGRAM 

Given  K  from  Part  I,  Part  II  program  performs  the  probe 
calculations  of  the  corresponding  initial  approach  angles 
(<Jj0)a>  ^o^B'  ^o^  *'  and  ^o^NL  based  on  the  relationship 


t(»o>  -  KA,0,.,NL[Sin(»o)A,B,*,NI,)P':1  . 

“ro^tl+(COS^o)A,Bi*,NL>  ^  =  0 

The  program  used  a  fixed  torpedo  speed  of  VT  =  35  kts  (18  m/s) 
and  an  initial  detection  range  of  rQ  *  R^  »  250  m.  The  purpose 
of  this  general  probe  program  is  to  locate  roots  of  the  given 
function  y  =  f(<t>Q)  to  evaluate  the  slope  of  the  tangent  line, 
and  to  find  the  maximum  and  minimum  points  on  a  graph.  We 
will  only  use  the  first  feature. 
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The  program  requires  a  subroutine  for  the  function  to 
be  investigated.  This  subroutine  starts  at  program  location 
140  and  is  located  at  label  2nd  C.  The  user's  instructions 
for  the  program  are  as  follows: 


Procedure 

Enter 

Press 

Display 

Enter  data 

Starting  Value 
of  ' 

2nd  D 

*o 

Enter  data 

The  increment  A <b 

o 

2nd  E 

A*0 

Calculate 

*<VMo> 

— 

A 

f<VA*o) 

Calculate 

f<V4V 

— 

B 

l-h 

-e- 

0 

1 

c> 

-e- 

0 

Display  current 
value  of  <J>0 

— 

E 

d> 

To 

If  the  value 

of  A$o  is  chosen  too 

large , 

A(J>q  may  be 

replaced  by  A$q/10  by  pressing  label  2nd  A.  Similarly,  if 
a  larger  value  of  A$o  is  required,  A<j>Q  can  be  replaced  by 
10A<j>o  by  pressing  label  2nd  B.  A  listing  of  the  program 
steps  follows,  giving  location  (LOC) ,  code  (COD),  key  symbol 
(KEY),  and  comments. 

Associated  with  the  subroutine,  it  should  be  noted  that 
the  value  of  K,rQ  and  p  are  entered  separately  in  the 
memory  locations  10,  11,  and  12,  respectively. 
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Fig.  2.  Location  of  Weather  Stations  Along  the  Norwegian 
Coast. 


Per  Year 


1  23456789 


Sea  State 


Fig.  4.  Histogram  of  Significant  Wave  in  Percent  per  Year 
from  the  Weather  Station  "Andenes." 
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j.  7.  Resonant  Bubble  Density  at  120  kHz  as  a  Function 
of  Depth. 

138 


A'***  t  r  *4  /Xfc* 


Sound  Pressure  Spectrum  Level  (dB  re  .0002  dyne/cm 


Pig.  8.  Ambient  Noise  Level  Curves. 
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Change  in  Sound  Speed  (m/s) 


Summer  4,5 
.Spring  6,7 


Range  (m) 


Frequencies  of  Wind  Directions  in  Percent  for 
Stations  Along  the  Coast  from  "Hillesoy  to 


Ona 


Absorption  Coefficient  (dB/m) 


Frequency  (kHz) 


Pig.  14.  Absorption  Coefficient  in  dB/ro  as  a  Function  of 
Frequency . 
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Fig.  17.  Spectrum  Level  as  a  Function  of  Total  Gas  Content 
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Fig.  20.  Required  Input  S/N  Ratio  and  BT  Product  for  Various 
Operating  Probabilities. 
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Fig.  21.  Probability  of  Detection  Versus  Input  S/N  Ratio  for 
Various  Thresholds. 
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eometry  at  the  Sea  Surface  Scattering. 
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Given  •  h4  ■  Source  depth 

hr  -  receiver  depth 
^<t>  *  haU  beam  width  of  receiver 
W  s  wind,  speed 
Varying:  R0  detection  range 


Geometry  : 
•  0 

•  a 

■  b 

•  R 


From:0[  From  .a  ,b 

R  A 


From:  0,<J‘a> 

•  a 

•  v 

•  firf(v} 

•  erfc(v)»l-erfv 

•  B 
•5(9) 


From:  W 


<  *5‘a> 


<&a>-  A  S(a>  Jt* 
<PP*>  R1l8n<TJ> 


I^Ovreci 

■“  lL0trect 


Fig.  26.  Calculation  Scheme  of  the  Surface  Scattering  Effect. 
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a) 

Fig.  29.  a) 
b) 


b) 

Chains  of  Successive  Scattering. 

Scattering  Pattern  Going  Through  the  Same 
Scatterer  More  than  Once. 
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Wave 


Fig.  30.  Scattering  Pattern. 
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Cross  Section  (m 


163 


n(a)/da  (m-  in  a  1  pm  band) 


Cross  Section  (m 


Pig.  32.  a  and  a  for  30  kHz 


n(a)/da  (m_  in  a  1  pm  band) 


Fig.  34.  a.  (a) n  (a)  da  for  60  kHz 
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crn(a)da  (m~  in  a  1  pm  band) 


f  *30  kHz 
hfaR')  *20 


Z  »  6  m 


Bubble  Radius  (pm) 


Fig.  37.  o„(a)n(a)da  for  30  kHz. 
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